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ABSTRACT 



The goal of this thesis is the development and implementation of a non-pcrturbative solution 
method for Wegner's flow equations. We show that a parameterization of the flowing Hamiltonian 
in terms of a scalar function allows the flow equation to be rewritten as a nonlinear partial 
differential equation. The implementation is non-perturbative in that the derivation of the PDE 
is based on an expansion controlled by the size of the system rather than the coupling constant. 
We apply this method to the Lipkin model and obtain very accurate results for the spectrum, 
expectation values and eigenstates for all values of the coupling and in the thermodynamic limit. 
New aspects of the phase structure, made apparent by this non-perturbative treatment, are 
also investigated. The Dicke model is treated using a two-step diagonalization procedure which 
illustrates how an effective Hamiltonian may be constructed and subsequently solved within this 
framework. 



OPSOMMING 

Die doel van hierdie tesis is die ontwikkeling en implementering van nie-steuringsteoretiese 
metodes vir die oplos van Wegner se vloeivergelykings. 0ns toon dat 'n parameterisasie van 
die vloeiende Hamilton operator in terme van 'n skalare funksie daartoe lei dat die vloeivergelyk- 
ing herskryf kan word as 'n nie-liniere parsiele differensiaal vergelyking. Hierdie implementering 
is nie op steuringsteorie gebaseer nie, en maak slegs gebruik van 'n reeks uitbreiding wat deur die 
grootte van die sisteem, eerder as die koppelingsterkte, beheer word. 0ns pas die metode toe op 
die Lipkin model en verkry akkurate resultate vir die spektrum, vcrwagtingswaardes en eietoes- 
tande vir alle waardes van die koppelingsterkte cn in die termodinamiese limiet. Nuwe aspcktc 
van die fase struktuur wat deur hierdie bcnadcring na vore gebring is, word ook bcsprcck. Die 
Dickc model word ondersoek deur 'n twee-stap diagonalisasie prosedure wat dien as voorbeeld 
van hoe effektiewe operatore binne hierdie raamwerk gekonstrueer en opgelos kan word. 
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Introduction 



The renormalization group and its associated flow equations U have become an indispensable 
tool in the study of modern physics. Its applications range from the construction of effective 
theories to the study of phase transitions and critical phenomenon. It also constitutes one of the 
few potentially non-perturbative techniques available for the treatment of interacting quantum 
systems. Our interest lies with a recent addition to this framework proposed by Wegner and 
separately by Glazek and Wilson 0, namely that of flow equations obtained from continuous 
unitary transformations. The flow equation in question describes the evolution of an operator, 
typically a Hamiltonian, under the application of a sequence of successive inflnitesimal unitary 
transformations. These transformations are constructed so as to steer the evolution, or flow, of 
the operator towards a simpler, possibly diagonal, form. The major advantage of this approach 
is that no prior knowledge of these transformations is needed as they are generated dynamically 
at each point during the flow. These attractive properties have led to applications to several 
diverse quantum mechanical problems, including that of electron-phonon coupling |2], boson and 
spin-boson models P| I42[ H3] . the Hubbard model the Sine-Gordon model and the Foldy- 
Wouthuysen transformation [Jj. The Lipkin model has also been particularly prominent among 
applications [HI El HOI IIH I12j . More recently the flow equations have been used to construct 
effective Hamiltonians which conserve the number of quasi-particles or elementary excitations 
in a system Examples of quasi-particles treated in this manner include triplet bonds on 

a dimerized spin chain jl4| I15| I17j and particle-hole excitations [1H\ I19j in Fermi systems. 
Studies of these equations in a purely mathematical context, where they are known as double 
bracket flows, have also been conducted j2()| I21j. 

The versatility of this approach should be clear from the references above. Unfortunately 
the practical implementation is hampered by the fact that the Hamiltonian typically does not 
preserve its form under the flow, and that additional operators, not present in the original 
Hamiltonian are generated |^ H] . In general we are confronted with an inflnite set of coupled 
nonlinear differential equations, the truncation of which is a highly non-trivial task. Perturbative 
approximations do allow one to make progress, although the validity of the results is usually 
limited to a single phase. 

It is the aim of this thesis to develop methods for the non-perturbative treatment of the 
flow equations. We will do so via two routes, both of which involve the representation of the 
flowing Hamiltonian as a scalar function and a systematic expansion in where N represents 
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the system size. This allows the original operator equation to be rewritten as a regular partial 
differential equation amenable to a numeric or, in some instances, analytic treatment. The bulk 
of this work comprises of the detailed application of these methods to two simple but non-trivial 
models. These calculations are found to reproduce known exact results to a very high accuracy. 

The material is organized as follows. Chapter 1 provides a brief overview of the flow equation 
formalism with particular emphasis on the role of the generator. In Chapter 2 we present the two 
solution methods. The first is based on an expansion in fluctuations controlled by the system 
size, while the second makes use of non-commutative coordinates to rewrite the flow equation 
in terms of the Moyal bracket [22] • The application of these methods to the Lipkin model 
constitutes the third chapter. We are able to calculate both eigenvalues and expectation values 
non-perturbatively and in the thermodynamic limit. New aspects of the phase structure, made 
apparent by the non-perturbative treatment, are also investigated. These results, published in 
[2^\ . have subsequently led to further studies of the phase structure in |21]. In Chapter 4 the 
Dicke model is treated using a novel two-step diagonalization procedure. Although we derive the 
flow equation non-perturbatively its complexity necessitates a partially perturbative solution. 
Despite this our approach serves as a valuable example of how an effective Hamiltonian may be 
constructed and subsequently solved within this framework. 



viii 



CHAPTER 1 
Flow equations 



1.1 Overview 

The central notion in Wegner's flow equations is the transformation of a Hamiltonian H 
through the application of a sequence of consecutive infinitesimal unitary transformations. It is 
the continuous evolution of H under these transformations that we refer to as the flow of the 
Hamiltonian. These transformations are constructed to bring about decoupling in leading 
to a final Hamiltonian with a diagonal, or block diagonal, form. The major advantage of this 
approach is that the relevant transformations are determined dynamically during the flow, and 
no a priori knowledge about them is required. 

We are led to consider a family of unitary transformations U{t) which is continuously parametrized 
by the flow parameter i S [0, cxd). U{1) constitutes the net effect of all the infinitesimal trans- 
formations applied up to the point in the flow labelled by L At the beginning of the flow f7(0) 
equals the identity operator. The evolution of U{1) is governed by 

^ = -u{i)v[i) (1.1) 

where rj{t) is the anti-hermitian generator of the transformation. Applying U{t) to H produces 
the transformed Hamiltonian H{t) = U\tjHU{t) for which the flow equation reads 

'^-^ = wim]- (1.2) 

To be consistent in the calculation of expectation values the same transformation needs to be 
applied to the relevant observables. An eigenstate of H transforms according to 

\cP,i)=U\£)\cP), (1.3) 

while the flow of a general observable 0{i) = {£)OU{£) is governed by 

^ = m,om. (1.4) 

Expectation values in the original and transformed basis are related in the usual way: 

(<A|O|0) = {cP\U{£)U\i)OUm\£M = (0,^|O(£)|(A,^) . (1.5) 

1 
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1.2 The role of T]{i) 

Much of the versatihty of the flow equation method stems from the freedom that exists in 
choosing the generator ri(i). Several different forms have been employed in the literature, and 
we will explore the consequences of some of these next. For the moment we restrict ourselves to 
the finite dimensional case. 

1.2.1 Wegner's choice 

In Wegner's original formulation '.2 r}{l) was chosen as the commutator of the diagonal part 
of H{t), in some basis, with H{t) itself, i.e. 

ll{t) = [Hdiag{i).H{i)]. (1.6) 

It was shown that in the £ — > oo limit H[t} converges to a final Hamiltonian H{co) for which 

[Hdiag{'X^).H{^)]=i). (1.7) 

We conclude that the effect of the flow is to decouple those states which correspond to differing 
diagonal matrix elements. In general this leads to a block-diagonal structure for H(oo). 

We will not use this formulation as other choices exist which offer greater control over both 
the type of decoupling present in H(oo) (i.e. the fixed point of the flow) and the form of H(i) 
during flow (i.e. the path followed to the fixed point). 

1.2.2 rj{e) = [G,Hii)] 

An alternative to Wegner's formulation is 

v{e) = [G,H{e)] (1.8) 

where G is a fixed (^-independent) hermitian operator of our choice. It is straightforward to 
show that H{i) converges to a final Hamiltonian which commutes with G. The proof rests on 
the observation that 

f^tT{H{e) - Gf = -2tr([G,F(£)]t[G,F(£)]) < 0, (1.9) 

where the positivity of the trace norm has been used. It follows that tr(ff(£) — G)^ is a monoton- 
ically decreasing function of i that is bounded from below by 0, and so its derivative must vanish 
in the ^ ^ oo limit. The right-hand side of (|1.9|) is simply the trace norm of r](i) = [G,H{i)], 
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and so we conclude that r/(oo) = 0. Choosing a diagonal G clearly leads to a block-diagonal 
structure for H{oo) where only states corresponding to equal diagonal matrix elements of G are 
connected. Put differently, G assigns weights to different subspaces through its diagonal matrix 
elements. The flow generated by r]{£) then decouples subspaces with differing weights. In par- 
ticular, a non-degenerate choice of G will lead to a complete diagonalisation of H. Furthermore, 
it can be shown [20] that the eigenvalues of H, as they appear on the diagonal of H{oo), will 
have the same ordering as the eigenvalues (diagonal matrix elements) of G. We can summarize 
this by saying that the flow equation generates a transformation that maps the eigenstates of H 
onto the eigenstates of G in an order preserving fashion. 

It is worth noting that this ordering can only take place within subspaces that are irreducible 
under G and H. The reason for this is that the flow equation clearly cannot mix subspaces that 
are not connected by either H or G. We will see several examples of this later on. 

1.2.3 The structure preserving generator 

Whereas the formulation above provides a good deal of control over the structure of H{oo), the 
form of H{i) at finite £ is generally unknown. This is due to non-zero off-diagonal matrix elements 
appearing at finite i that are not present in either the initial or final Hamiltonian. For example, 
a band-diagonal Hamiltonian may become dense^ during flow and still converge to a diagonal 
form. From a computational point of view it is clearly desirable that H{i) assumes as simple 
a form as possible. Of particular interest in this regard are generators which preserve a band 
diagonal, or more generally band block diagonal structure present in the original Hamiltonian. 
These types of generators have been applied to a wide range of models El El El El E] > and 
are particularly attractive in that they allow for a clear physical interpretation of the transformed 
Hamiltonian H{oo). 

First we introduce an operator Q with integer eigenvalues which will serve as a labelling device 
for different subspaces in the Hilbert space. Associated with each distinct eigenvalue qi of Q is 
the corresponding subspace of eigenstates Qi where i = 1,2,3,.... We also assume, without 
loss of generality, that qi < gj+i for all i. The Hamiltonian H is said to possess a band block 
diagonal structure with respect to Q if there exists an integer such that {i \ H\ j) = for all 
\i) G Qi and \j) G Qj whenever \qi — qj\ > N. This selection rule clearly places a bound on 
the amount by which H can change Q. The matrix representation of such an operator in the Q 
basis typically has a form similar to that shown in Figure ITTT] fa). For the cases we will consider, 
dense matrix is one of which the majority of elements is non-zero. 
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Figure 1.1: The typical structure of (a) band block diagonal and (b) block diago- 
nal matrices. Darker shaded blocks correspond to matrix elements which connect 
increasingly distant Q-sectors. Unshaded areas are filled with zero elements. For 
purposes of illustration we used = 2 in (a). 



and for those treated in the literature, it is possible to group terms in the Hamiltonian together 
based on whether they increase, decrease or leave unchanged the value of Q. This leads to the 
form 

n=N n=-N 

H = To+Y,Tn+ Yl (1-10) 

n=l n=— 1 

T+ T- 

where Tn changes Q by n, i.e. [Q,T„] = nT„. Clearly Tq is hermitian, while T„ and T_„ are 
conjugates. Tq is responsible for scattering within each Q-sector. 

Flow equations are used to bring this Hamihonian into a form which conserves Q, i.e. for 
which [H(oo),Q] = 0. This effective Hamihonian wih be block diagonal, similar to Figure ITTT] 
(b), with each block containing new interactions generated during the flow. We require that the 
band block diagonal structure of H is retained at finite £, and so the flowing Hamiltonian should 
be of the form 

n=N n=-N 

H{e) = To{£) + ^ Tn{e) + Tn{l) . (1.11) 

n=l n=— 1 

T+(£) T_{<?) 

The generator which achieves this is 



r/(£) = r+(£)-r_(^). 



(1.12) 
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The corresponding flow equation reads 

dH{e) 



2[r+w,r_(^)] + [T+(^) -r_(£),ro(^)] . (1.13) 



The combinations of T+,T_ and Tq appearing in (|1.13|) clearly cannot generate scattering be- 
tween Q-sectors differing by more than A^, and so the band block diagonal structure of H{Q) is 
preserved. It can be shown that this generator guarantees convergence to the desired form 
of H{oo) which conserves Q. 

Before proceeding, let us point out why the choice of the previous section r]{i) = [Q,H{£)] 



does not lead to a flow of the form although it does produce the correct flxed point. Using 

the form of H{0) we see that 

r,{0) = [Q,H]=T',-TL (1.14) 



where Tj_ = Yl^=i '^'^±n-, and so 

dH{(.) 



M 



= [r;-r:,ro + r+ + r_] . (1.15) 

1=0 



Whereas commutators between terms which increase (decrease) Q dropped out in (|1.13|) this is 
not the case in (|1.15|) . In general [7±,T-j-] will contain terms which change Q by up to 2N — 1, 
thus destroying the band block diagonal structure. This process will continue until, at flnite ^, 
H{t) will connect all possible Q-sectors. 

Although our treatment has dealt with Q largely in abstract terms, this is not generally the 
case in the literature. Previous applications of this method present a physical picture of Q 
as a counting operator for the fundamental excitations describing the low energy physics of a 
system. Examples of these include triplet bonds on a dimerized spin chain jl4| II 5| I16L Wf\ and 
quasi-particles (particle-hole excitations) |18| llHj in Fermi systems. The flow equation is used to 
obtain an effective Hamiltonian that conserves the number of these excitations. These "eff'ective 
particle-conserving models" ^3] have been studied in considerable detail within a perturbative 
framework. 

Finally we point out that the approaches of Sections II . 2.21 and II . 2 . 31 coincide when H only con- 
nects sectors for which Q differs by some fixed amount N , i.e. {i \H\ j) = when \qi — qj\ ^ {0, A^}. 
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1.3 Flow equations in infinite dimensions 

Proofs concerning the convergence properties of the flow equations discussed thus far rely 
on the invariance of the trace under unitary transformations. In infinite dimensions the trace 
is not generahy well defined, and the question of convergence depends on special properties 
of the Hamiltonian. The most important property in this regard is the boundedness of the 
spectrum. It has been shown that all the methods described above will converge provided that 
the Hamiltonian possesses a spectrum bounded from below [21 E]- Aspects of the flow equations 
in infinite dimensions were investigated in |25j and |2S] using the bosonic Hamiltonian 

H = a^a + ^(^a^ + a^'^Y (1.16) 

and we digress for a moment to discuss this case in more depth. It is well known that for 
a < 1 this model possesses a harmonic spectrum which can be found by applying the Bogoliubov 
transformation. When a > 1 the spectrum forms a continuum which is unbounded from both 
above and below. We will contrast the behaviour of the flow equations resulting from Wegner's 
choice of ri{£) = [Hdiag{(), H{£)] and the choice r/(£) = [o^^o, H{£)'j . In both these cases H{£) has 
the simple form 

Hii) = a{e) + u;{e)a^a + ^ (a^ + a^') , (1.17) 

and the flow equation closes on a coupled set of three differential equations for cr{i), uj{i) and 
a{i). When < a < 1 in the starting Hamiltonian we find that both choices of rj lead to the 
fixed point 

H{oo) = ^{uj{oo) - a) +uj{oc)a^a (1.18) 

where w(oo) = \/l — a'^; a result matching that of the Bogoliubov transformation. 

Next, consider the a > 1 case. When we solve the flow equation for r]{£) = [a^a, H{£)~^ 
we no longer find convergence. Surprisingly, Wegner's choice for the generator still produces a 
convergent flow, although not one leading to a diagonal form: 

H{oo) = -l + \\/o? - 1 (a^ + at') . (1.19) 

Note that this does not contradict the assertion made earlier that [i^^jagCoo), //(cxd)] = 0. 

Admittedly this example does represent an extreme case. For physically sensible Hamiltonians 
the spectrum is always bounded from below, as one would expect in order to have a well-defined 
ground state. For the models that we will consider the flow equations may be safely applied 
without further modification. 



CHAPTER 2 
Solving the flow equation 



We have seen the theoretical capabihties of flow equations with regard to the diagonalization 
of Hamiltonians and the construction of effective operators. The practical implementation of 
this method is, however, hampered by the difficulty of solving the resulting operator differential 
equation. On the operator level this is due to the generation of additional operators during the 
flow that were not present in the original Hamiltonian. This leads to an extremely large set of 
coupled differential equations for the coupling constants of these terms. Generally some kind of 
approximation is required in order to continue. The usual approach consists of replacing H{i) 
by a simpler parametrized form for which the flow equation closes on a set of equations of a 
tractable size. A particular parametrization is usually selected on the basis of a perturbative 
approximation, or by using some knowledge of the relevant degrees of freedom in the problem. 
In general this approach is only valid for a limited range of the coupling constant, and tends to 
break down when the system exhibits non-perturbative features, i.e. non-analytic behaviour in 
the coupling constant. 

We will introduce two new approaches to this problem which allow us to treat the flow equation 
in a non-perturbative way. The first involves a systematic expansion in fluctuations, which is 
controlled by the size of the system rather than the coupling constant. The second makes use of 
non-commutative variables to recast the flow equation as a regular partial differential equation. 
This approach also involves an expansion controlled by the system size. Although there are 
similarities between these two methods we will treat them separately and then show how they 
produce the same results in specific cases. 

2.1 Expansion in fluctuations 

We consider the flow of a time reversal invariant hermitian operator H(i) which acts on a 

Hilbert space 7i. We assume 7i to be finite dimensional since the expansion ()2.1() below is easily 

proved in this case (see Appendix^. This restriction is, however, by no means essential. If 

the Hilbert space is infinite dimensional the expansion (|2.1|) still applies to bounded operators 

j27j and subsequently one may approach the infinite dimensional case by studying the flow of 

a bounded function of the Hamiltonian, instead of the Hamiltonian itself. Alternatively one 

can introduce a cut-off in some basis, e.g., a momentum cut-off and study the behaviour of the 

system as a function of the cut-off. 

The aim of this section is to develop a parameterization of the flowing Hamiltonian which 

7 
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allows for a systematic expansion controlled by the fluctuations, rather than the coupling con- 
stant. This is in the same spirit as the semi-classical expansion of quantum mechanics, based 
on an expansion in orders of h, and corresponds to resumming certain classes of diagrams in the 
perturbation series of the coupling constant. 

Let Hq and Hi be hermitian operators acting on Ti. which together form an irreducible set. 
What follows holds for any irreducible set, however, the case of two operators appears naturally 
in flow equations as the Hamiltonian can usually be written in the form H = Hq + Hi where 
the spectrum and eigenstates of Hq are known. Note that if Hq and Hi are reducible on TC, but 
irreducible on a proper subspace of 7i, the problem can be restricted to this smaller subspace 
making the irreducibility of Hq and Hi a very natural requirement. In Appendix ^ it is shown 
that any operator acting on 7i can be written as a polynomial in Hq and Hi. In particular this 
holds for H(i), and so we may write 

Hii) = r + VHi + r^HiHj + V^^HiHjHk + ... (2.1) 

where each r*-'"' coefficient is a function of £ and repeated indices indicate sums over and 1. 
Since H(i) is both hermitian and real it follows that the F's are invariant under the reversal of 
indices: 

r"i«2...nfc ^ pnfenfc_i...ni > Q, G {0, 1}. (2.2) 

Now define AHi = Hi — {Hi) as the fluctuation of Hi around the expectation value {Hi) = 
{(l)\Hi \ (f)), where {(p) is some arbitrary state. By setting Hi = AHi + {Hi) in equation (|2.1j) we 
obtain H{i) as an expansion in these fluctuations 

H{i) = f + VAHi + r^AHiAHj + V^^AHiAHjAHk + ..., (2.3) 

where 

f = T + r{Hi)+r^ {H>){Hj)+r^^{H,){H,){Hk) + ... (2.4) 

r = r + (r^' + f^) {h^) + {v^^^ + + r'^^'*) {Hj) {Hk) + ... (2.5) 
= + (r'''^ + r''^^' + r^^'^) {Hk) + ... (2.6) 

and so forth. Note that f, when viewed as a function of i, {Hq) and {Hi), encodes information 
about all the expansion coefficients F*-' " appearing in equation l\2.1^ . A natural strategy that 
presents itself is to set up an equation for f as a function of £, {Hq) and {Hi) in some domain V 
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of the {Ho)~{Hi) plane, determined by the properties of the operators Hq and Hi. In particular 
we require that this domain includes values of (^^o) ranging from the largest to the smallest 
eigenvalues of Hq. For this purpose we require a family of states \(j){X)), parameterized by a 
continuous set of variables A E P, such that as A is varied over the domain 1), (Hq) and (Hi) 
range continuously over the domain D. A very general set of states that meets these requirements 
are coherent states [2Hj. In this representation we may consider {Hq) and {Hi) as continuous 
variables with each T a function of {Hq), {Hi) and the following relations hold generally: 

and = -— ilA^. (2.7) 



d{Hi) 2d{H,)d{Hj)' 

Coefficients with three or more indices cannot be written in this way, since, for example, T^^^ 
need not be equal to f^^. The general relationship between the true coefficients and derivatives 
of r are 

an+mp 

n!m! ^ fP^"'™) (2.8) 



d{Ho)-d{Hi)^ 



where the sum over p is over all the distinct ways of ordering n zeros and m ones. For example 

1 d^r 1 



31 d {Ho)' d {Hi)' 3 



(pOOi ^ f 010 ^ f 100) _ (2 9) 



Replacing T^^' by the left-hand side of the equation above is equivalent to approximating F''^^ 
by the average of f ooi^f and F^*''^. In general this amounts to approximating f by the 
average of the coefficients corresponding to all the distinct reorderings of ni, . . . , n^. 

To set up an equation for T{{Hq) , {Hi) , £) we insert the expansion (|2.:-{j) into the flow equation 
with r](i) = [Hq, H(i)] and take the expectation value with respect to the state |(/)(A)). The left- 
and right-hand sides become a systematic expansion in orders of the fluctuations {AH^ ■ ■ ■ AHn^ ) ■ 



IJL (71 — ■ — • 

— + ^{AHAH,) + ... = VV{[[AHq,AH,\,AH,]) 

+ fv^^^ {[[ahq, ah,], ah ah J,]) 

+ f''^P{[[AHQ,AH,,AHi,],AHj]) (2.10) 

+ fiii2fhn ([[^^^^ AHi,AHQ,AHj,AH,,]) 



Note that writing Hq or AHq in the ffist position of the double commutator on the right-hand side 
of equation (|2.1())) is a matter of taste. The expectation values appearing above are naturally 
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functions of A and may be written as functions of (Hq) and {Hi) by inverting the equations 
(Hi) = (0(A) \Hi\ (/>(A)) to obtain X{{Hq) , {Hi)). With an appropriately chosen state, such as a 
coherent state, the higher orders in the fluctuation can often be neglected, as the expansion is 
controlled by the inverse of the number of degrees of freedom (see Appendix O for an explicit 
example). A useful analogy is the minimal uncertainty states in quantum mechanics which 
minimizes the fluctuations in position and momentum. The choice of state above aims at the 
same goal for Hq and Hi. Clearly it is difficult to give a general algorithm for the construction 
of these states and this property has to be checked on a case by case basis. If this is found to 
be the case we note that replacing the r^-^ ^'s with more than two indices by a derivative will 
introduce corrections in (|2.1()|) of an order higher than the terms already listed, or, on the level 
of the operator expansion, corrections higher than second order in the fluctuations. Working 
to second order in the fluctuations we can therefore safely replace T^^--- by the derivatives of a 
function / {{Hq) , {Hi) ,£) = r{{Ho) , {Hi),e) and write for H{e): 

Hii) = f + f'AHo + f'AHi + ^fH^Ho,AHi} + ^f^AH^ + ^f'AHf + ... 

= ^ — — ■ ^ (Distinct orderings of i AHq's and j AH[s^ (2-11) 

i,j=o (^ + -^')- 

where {•,•} denotes the anti-commutator and 

^ ^_ .2.12) 

d{Ho)'d{Hi)^ 

This turns the flow equation 1)2.10(1 into a nonlinear partial differential equation for / {{Hq) , {Hi) , i), 
correct up to the order shown in (|2.10|) . The choice of coherent state, the corresponding calcu- 
lation of the fluctuations appearing in equation (|2.10j) and the identification of the parameter 
controlling the expansion are problem specific. There are, however, a number of general state- 
ments that can be made about the flow equation and the behaviour of / {{Hq) , {Hi) The 
first property to be noted is that since H{oo) is diagonal in the eigenbasis of Hq it should become 
a function of only Hq, provided that the spectrum of Hq is non-degenerate. This is reflected in 
the behaviour of / by the fact that / {{Hq) , {Hi) , £ = oo) should be a function of {Hq) only. This 
is indeed borne out to high accuracy in our later numerical investigations. This function in turn 
provides us with the functional dependence of H{oo) on Hq. Keeping in mind the unitary con- 
nection between H(oo) and H this enables us to compute the eigenvalues of H straightforwardly 
by inserting the supposedly known eigenvalues of Hq into the function / {{Hq) , cxd). 

A second point to note is that the considerations above apply to the flow of an arbitrary her- 
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mitian operator with time reversal symmetry. It is easily verified that the transformed operator 
0{i) = U\£)O{0)U{i) satisfies the flow equation: 

^ = m,0{£)] (2.13) 

where the same choice of r]{i) as for the Hamiltonian has to be made, i.e., r]{i) = [Hq, H{i)]. The 
expansion (|2.11() can be made for both operators H(l) and 0{£). Denoting the corresponding 
function for 0{i) by g {{Hq) , (-f^i) , i) the flow equation (|2.13j) turns into a linear partial differen- 
tial equation for g {{Hq) , (Hi) , £) containing the function / ((^^o) ) {Hi) , £), which is determined 
by (|2.10|) . The expectation value of O = 0(0) with respect to an eigenstate \En) of H can be 
expressed as 

{En\0\En) = (En U{£)U^{£)OU{£)U\£) E^) = {En,£\0{£)\ E^J) (2.14) 



where \En,£) = U^{£) \En)- In the limit £ ^ cc the states |ii^„,oo) are simply the eigenstates 
of Hq, which are supposedly known. In this way the computation of the expectation value 
{En \0\ En) can be translated into the calculation of expectation values of the operator 0{oo), 
obtained by solving the flow equation (|2.13)) . in the known eigenstates of Hq. 



2.2 Moyal bracket approach 

Next we present an approach based on the Moyal bracket formalism |22| . Let Ti. denote the 
D-dimensional Hilbert space of the Hamiltonian under consideration. We define two unitary 
operators h and g that act irreducibly on 7i and satisfy the exchange relation 

hg = e-'^gh. (2.15) 

Since g is unitary its eigenvalues are simply phases. Let e**^ be one such eigenvalue, and consider 
the action of gh on the corresponding eigenstate \a): 

gh \a) = e'^hg \a) = e''-^+''^h \a) . (2.16) 

We see that h\a) is again an eigenstate of g with eigenvalue e*(^+'^). Since g and h act irreducibly 
on TC all the eigenstates of g can be obtained by the repeated application of to \(t). Furthermore, 
we may scale g so that it has an eigenvalue equal to one. 
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It follows that the eigenvalues and eigenstates of g take the form 

g\n) = e''^''\n) where n = 0, 1, 2, . . . , L> - 1 (2.17) 

while h acts as a ladder operator between these states: 

gh\n) = e''^hg\n) =e''^^''+^'^h\n) =^ h\n)(x\n + l). (2.18) 

The allowed values of 9 are found by taking the trace on both sides of h~^gh = e^^g, which leads 
to the requirement 

D-l 

iv{g) = e"' = 0. (2.19) 

n=0 

This fixes 9 at an integer multiple of 2tt/D. We choose 9 = 2tt/D as this ensures that g is 
non-degenerate, which is crucial for the construction that follows. 

In similar fashion to the previous section, we wish to represent flowing operators in terms of 
g and h. The main result in this regard is that the set 



V 



^D-^g^-h"" : n,m = 0,1,2, D - 1^ (2.20) 



forms an orthogonal basis for the space of linear operators acting on TC. The orthogonality of V 
follows from applying the trace inner product in the g'-basis to two members of V: 



D-l 

, „ ,4- „/ /1 _ r /„/ „\ 

tr 



(5^/1"*)^ g^'h^\ = (5^,^'tr [5("'-)J = 5^,^, e^''^"'-"^' = ^m^m'^u'D, (2.21) 

where bm,n equals one if m = n mod D and is zero otherwise. This, together with the observation 
that the dimension of the linear operator space equals \P\ = D"^, proves the claim. 
Consider two arbitrary operators U and V expressed in the V basis as 

U = Y Cn,mg''h^ 8.nd V=J2 C'n',m'9'''h"'' (2.22) 

n,m n',m' 

where Cn,m and C^, ^, are scalar coefficients. We use the convention of always writing the /I's 
to the right of the ^'s. The product of U and V then gives 

UV= C'„,^C;,,^,5'"+"*'/i"+"'e-^""^'^. (2.23) 



n.m.n' ,m' 
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Note the similarity in form between this product and the product of functions of regular com- 
muting variables. Only the phase factor, the result of imposing our ordering convention on the 
product, distinguishes the two. In fact, we may treat g and h as regular scalar variables pro- 
vided that we modify the product rule to incorporate this phase. Convenient variables for this 
procedure are a and /?, which are related to g and h (now treated as scalars) through g = e*" 
and h = e^^ . Having replaced operators by functions of a and (3 the modified product rule reads 

U{a, 13) * V{a, /?) = U{a, (3)e''^^0 ^"^(a, (2.24) 

where the a and (3 derivatives act to the right and left respectively. This is seen to be of the 
required form by using the fact that both g and h are eigenfunctions of da and dp: 

U{a,l3)e'^^^^"V{a,l3) = ^ C„,„ C7;,_„, e^'"" e*"^ efj^ e^™'° e^"'^, (2.25) 

n,m,n',m' g-iem'n 

which agrees with ()2.23|) . 

The *-operation is known as the Moyal product j22j . while the corresponding commutator 
[[/, y]^ = U*V — V*Uis the Moyal bracket. When H(i) and r]{£) are represented in this 
manner the flow equation becomes a partial differential equation in a, (3 and i: 

^Ei^llH = [^(«, (3, £), Hia, P, i)] , . (2.26) 

In its exact form this formulation is not of much practical value, since the operator exponent 
involved in the Moyal product is very difficult to implement numerically. A significant simpli- 
fication is achieved by expanding the operator exponent to first order in 9, which is known to 
scale like one over the dimension D of the Hilbert space. We expect this to be a very good 
approximation provided that the derivatives do not bring about factors of the order of D. This 
translates into a smoothness condition: we require that the derivatives of the relevant functions 
remain bounded in the thermodynamic limit as D goes to infinity. 

Using this approximation the Moyal product becomes, to leading order, 

U*V = UV + ieUpVa (2.27) 

while the Moyal bracket reads 



[U, VI = i9 {Up Va - Vp Ua) 



(2.28) 
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Partial derivatives are indicated by the subscript shorthand. The form of the flow equation is now 
largely fixed, up to the specific choice of the generator. As an example, consider the generator 
r]{a,l3) = [a,H(a, (3)]^ which we will use later on. In this case the flow equation becomes 

^ = 9^ {HppHa — HfsHfSa) ■ (2.29) 

The remaining problem is that of constructing the initial conditions, i.e. H{0), in terms of 
g and h (or equivalently a and /?) in such a way that the smoothness conditions are satisfied. 
The reader may have noticed that we have not specified how the realization of g and h should 
be constructed on TC. Put differently, there is no obvious rule which associates a specific basis 
of TC with the eigenstates of g. It seems reasonable that this freedom may allow us to construct 
smooth initial conditions through an appropriate choice of basis, whereas a malicious choice 
could produce very poorly behaved functions. We know of no way to proceed on such general 
terms, and we will instead tackle this problem on a case-by-case basis. In all of these we will 
use the algebraic properties of operators appearing in the Hamiltonian to reduce this problem 
to one of representation theory. 



Although two operators are clearly the minimum required to construct a complete operator 
basis, it is also possible to introduce multiple such pairs. This would be a natural choice when 
TC is a tensor product of Hilbert spaces TCi (i = 1, . . . , m), each of which is of a high dimension. 
We can introduce m pairs of operators {gi, hi} which satisfy higi = e~^^^gihi and higj = gjhi for 
all i ^ j- In the same way as before this leads to m pairs of scalar variables {ai,Pi} for which 
the product rule, to first order in the 0j's, is 

m 

U{cx,f3)*V{a,P) = UV + iY,OjUp,Va^. (2.30) 

Note that {ai, (3i} are analogous to conjugate position and momenta coordinates representing the 
independent degrees of freedom of the system. The Moyal bracket acts like the Poisson bracket 
for these coordinates: 



[ai,Pj]^ = -i65ij and [aj,aj]^ = = 0. (2.31) 

This formulation strongly suggests an analogy with semi-classical approximation schemes. Our 
approach to solving the fiow equation is indeed very closely related to the Wigner-Weyl-Moyal 
[221 129j formalism, which describes the construction of a mapping between quantum operators 
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and functions of classical phase space coordinates. This allows for the description of a quantum 
system in a form formally analogous to classical dynamics. When applied to the flow equations 
this formalism produces results similar to those obtained before. The central approximation 
again involves a non-perturbative expansion, but which is now controlled by h, and so is semi- 
classical in nature. Let us formalize some of these notions in the context of a single particle 
in 3 dimensions. The relevant Hilbert space is 7i = L'^(R^) and the position and momentum 
operators satisfy the standard commutation relations 



[xi,pj] = ihSij i,j = 1,2,3. 



(2.32) 



We first introduce the characteristic operator [201 



U{t, s) = e'^ Pe'^-^ 



(2.33) 



where t = (ti,t2;^3) £ x = {xi,X2,x^) and similar for s G and p. Varying the arguments 
of C/(t, s) over their domains produces a set of operators analogous to V (equation (|2.2flj) l. where 
the discrete powers n and m correspond to the continuous labels t and s. We again find both 
completeness and orthogonality with respect to the trace norm: 



tr 



U{t',s')^U{t,s) 



dx 



X 



^ [ da;dpe*P (*-*')e*^ (^-^') 
rhy J 



{2TTh) 



(2.34) 



Using this, an operator A(x,p) can be represented as 



A{x,p) = j Aids A{t,s)U{t,s) 



(2.35) 



where 



i(t, s) = tf' tr \u{t, s)'^A{x,p 



(2.36) 



is a scalar function. Now consider the product of two operators represented in this manner: 



A{x,p)B{x,p) = J dt ds dt' ds'A{t, s)B{t' , s')U{t + t' ,s + s')e 



(2.37) 



The non-commutativity of x and p gives rise to the scalar factor e which is the only 

element distinguishing this product from one of regular scalar functions. We conclude, as before. 
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that the position and momentum operators may be treated as scalar variables provided that we 
modify the product rule to incorporate this phase. This leads to the Moyal product 

A{x,p) * B{x,p) = A{x,p)e''^^^-^^B{x,p) (2.38) 

where x,p G M"^ and dx • dp= dx^ dp.. Note that an expansion of the exponential is now 
controlled by h instead of 9 (x dim(7i)~^. To leading order the Moyal bracket is given by 

3 

[A{x,p),B{x,p)l = ihJ2 {Axi Bp, - Ap, B^,) , (2.39) 

1=1 

where the subscripts denote partial derivatives. 

We conclude that in a semi-classical approximation the Hamiltonian and generator may be 
replaced by scalar functions H{x,p, £) and r]{x,p, £), and that the flow equation is given in terms 
of the Moyal bracket by 

= [vix,p,e),H{x,p,£)l. (2.40) 



When solved to leading order in h this equation describes the renormalization of the Hamiltonian 
within a semi-classical approximation. Further quantum corrections can be included by simply 
expanding the Moyal bracket to higher orders in h. 



CHAPTER 3 
The Lipkin Model 



3.1 Introduction 

Since its introduction in 1965 as a toy model for two shell nuclear interactions the Lipkin- 
Meshov-Glick model j31j has served as a testing ground for new techniques in many-body physics. 
Here we will use it to illustrate both the working of the flow equations and the solution methods 
presented in the previous chapter. While the simple structure of the model will allow many 
calculations to be performed exactly, its non-trivial phase structure will provide a true test for 
our non-perturbative approach. 

We begin with an overview of the model, its features and the quantities we are interested in 
calculating. After pointing out some specific aspects of the flow equations for the Lipkin model 
we proceed to treat the equations using the methods developed earlier. Finally we present the 
results obtained from the numerical solutions of the resulting PDE's and compare them with 
some known results. New aspects of the model, brought to the fore by this treatment, will also 
be discussed. 

3.2 The model 

The Lipkin model describes fermions distributed over two A-fold degenerate levels separated 
by an energy of For simplicity we shall take A^ to be even. Fermi statistics require that A^ < 2A, 
and accordingly the thermodynamic limit should be understood as A ^ oo followed by A^ — > oo. 
The interaction V introduces scattering of pairs between levels. Labelling the two levels by 
cr = ±1, the Hamiltonian reads 

^ = ^Y1 o-^.^^P.'^ ^ "p.'t^'.'tV -^Op (3-1) 

cr,p P,p' ,o- 

where the indices p and p' run over the level degeneracy 1 ... A. A spin representation for H can 
be found by introducing the su(2) generators 

Jz = ]^^(yal,,ap^a and J± = X] ^4,±i"p.t1- (3-2) 

a,p p 
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Together with the second order Casimir operator j"^ = j'^ + Jz + J-J+, these satisfy the regular 
su(2) commutation relations: 

[J„ J±] = ±J±, [J+, J_] = 2 J, and [J2, J±] = [J2, J,] = 0. (3.3) 

We divide H hy and define the dimensionless coupling constant A = NV/^^ to obtain 

H = J, + ^{Jl + Jl), (3.4) 

where all energies are now expressed in units of ^. The factor of in the definition of A 
brings about the 1/A^ in front of the second term, which ensures that the Hamiltonian as a 
whole is extensive and scales like N. Since ^H, J^] = the Hamiltonian acts within irreducible 
representations of su(2) where states are labelled by the eigenvalues of and Jz, i.e., \ j, m) = 
j{j + 1) \j, m) and Jz |i, m) = m \ j, m) for m = —j, . . . ,j. The Hamiltonian thus assumes a block 
diagonal structure of sizes 2j + 1. The low-lying states occur in the multiplet j = N/2, and 
we fix j at this value throughout, using the shorthand \m) = \j,m) for the basis states. The 
Hamiltonian can be reduced further by noting that it leaves the subspaces of states with either 
odd or even spin projection invariant. States belonging to one of these subspaces are referred 
to as having either odd or even parity. We denote the eigenstate of H with energy En by \En), 
where Eq < Ei < E2 ■ ■ ■■ When A = the ground state is simply |ii^o) = o\ _i . . . ajy _i |0) which 
is written as \Eq) = \ —j) in the spin basis. Non-zero values of A cause particle- hole excitations 
across the gap, and at A = ±1 the model exhibits a phase transition from an undeformed first 
phase to a deformed second phase. To distinguish the two phases we use the order parameter 

= 1 + (Jz) /j where (Jz) is the expectation value of Jz in the ground state. As we will show, 
is non-zero only within the second phase. 

The phases can be characterized further by the energy gap A = Ei — Eq which is positive 
in the first phase and vanishes like Vl — A^ as A approaches 1. In the second phase the ground 
states of the odd and even parity subspaces become degenerate, causing the parity symmetry to 
be broken and the corresponding energy gap to vanish. Further discussion of this model and its 
features can be found in pT] . 

3.3 The flow equation approach 

We will follow the formulations of Sections 11.2.21 and 11.2.31 and choose the generator as 



7^i£) = [Jz,H{e)] 



(3.5) 
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Let us consider the consequences of this choice. Firstly, since Jz is non-degenerate, we expect the 
final Hamiltonian H{oo) to be diagonal in the known spin basis. Furthermore the eigenvalues 
of H appear on the diagonal of H{oo) in the same order as in Jz, i.e. increasing from top to 
bottom. Secondly, note that H possesses a band diagonal structure with respect to Jz, which 
plays the role of Q in Section II. 2. 31 Based on our earlier discussion we expect this structure 
to be conserved during flow, and so H(£) will only connect states of which the spin projection 
differs by two. 

Let \En,() = U^{£) \En) denote the transformed eigenstates of H. From the ordering property 
of eigenvalues in H{oo) we conclude that \En, oo) = \—j + n), and so a general expectation value 
may be calculated through 

{En \0\ En) = {En, OO |0(oo)| En, oo) = {-j + u |0(oo)| - j + n) . (3.6) 

Before continuing we mention some of the previous applications of flow equations to the 
Lipkin model. The first of these was by Pirner and Friman 8 , who dealt with newly generated 
terms by linearizing them around their ground state expectation values. This yielded good 
results in the first phase but lead to divergences in the second; a common ailment of these types 
of approximations. Subsequent work by Mielke JHl relied on an ansatz for the form of H(i)^s 
matrix elements, while Stein jH] followed a bosonization approach. Dusuel and Vidal |lH ll2j used 
flow equations together with the Holstein-Primakoff boson representation to compute finite-size 
scaling exponents for a number of physical quantities. Their approach was also based on a 1/A^ 
expansion. A method which produced reliable results in both phases was that of j32j . This 
method also relied on the linearization of newly generated operators but, in contrast with [S], 
did so around a dynamically changing ("running") expectation value taken with respect to the 
flowing ground state. This "running" expectation value could be solved for in a self consistent 
manner, and then used in the flow equation for H{i). This lead to non-perturbative results for 
the ground state, excitation gap and order parameter in both phases. Although it described the 
low energy physics well, this method failed to produce the correct spectrum for the higher energy 
states. 

3.4 Expansion method 

In order to apply the method of Section r2.1l we require an irreducible set of operators {Hq, Hi} 
in terms of which the flowing Hamiltonian will be constructed. For the purposes of the Lipkin 
model we will choose Hq = Jz and ^^i = ^ ( J+ + J^) , i.e. the diagonal and off-diagonal parts 
of the Hamiltonian respectively. The reader may well remark that this does not constitute an 
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Figure 3.1: The domain of f{{Ho) , {Hi) , i). 



irreducible set on the entire Hilbert space, since the subspaces of odd and even spin projection 
are left invariant, as pointed out in Section [3.21 Indeed, it is only within a subspace of definite 
parity (i.e. odd or even spin projection) that this set is irreducible. However, we note that 
the flow equation will never mix these odd and even subspaces, and so the flow will proceed 
independently within each subspace. This ensures that our representation of H{£) will never 
require an operator which connects odd and even states, and it turns out that this choice of 
{i^Oi-f^i} is indeed sufficient. The same conclusion extends to U{£) by equation (|I.I|) . and then 
to any flowing operator 0{£) = {£)OU{£), provided that O can be represented in this way. 
We will calculate the averages of Hq and Hi with respect to the coherent state |28( 133) 



where z = rexp(i^) G C. See Appendix El for details on the properties of these states and the 
method by which the averages are calculated. We find that: 



This constitutes a mapping of the complex plane onto the domain shown in Figure 13.11 With 
these definitions in place we conclude that the flowing Lipkin Hamiltonian may be written as 



z) = (l+r2)-Jexp(zJ+)|-j) 



(3.7) 




(2j - l)r2 cos(26l) 
(l+r2)2 



(3.8) 



H{e) = f + f 



AHo + f'AHi + {Ai^o, AFi} + ^f^^AH^ + ^f^AHf + ■■■ (3.9) 



with Hq and Hi respectively the diagonal and off-diagonal parts of the original Hamiltonian. By 
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definition AHi = Hi — {z \ Hi \ z) and so / is defined on tlie domain pictured in Figure l3?Tl The 
initial condition becomes 

f{{Ho) , {Hi) ,£ = 0) = {Ho)+X {Hi) . (3.10) 

Due to the coherent nature of the state 1)3 .7(1 one might expect that terms corresponding to high 
order fluctuations in H3.9() will contribute less signiflcantly to {z\H{i)\ z) than the scalar term 
/. This statement can be made precise as follows: since the flowing Hamiltonian is extensive / 
should be proportional to j, as is the case with {Hq) and {Hi). To keep track of the orders of 
j we introduce the scaleless variables x = {Hq) /j, y = {Hi) /j and / = f /j. When taking the 
inner-product with respect to \z) on both sides of (|3.9j) the linear terms fall away and we obtain 

{Hil)) = f + {{AHo, AHi}) + (AHi) + (AHf) + • • • . (3.11) 

Each term is, up to a constant factor, of the form 

(Prod(i + k)) = (Prod(i + k)) (3.12) 

where Prod(« + k) denotes some arbitrary product of i + A; fluctuations. Using the results of 
AppendixOwe see that such a term is at most of order C)(ji~(*+'^)jl-(*+'^)/2J ). The leading order 
term corresponds to i + k = 0, i.e. the scalar term /. We conclude that / is the leading order 
contribution to {z \ H{i) \ z), expressed not as a function of z and z* , but rather of the averages 
{Ho) and {Hi). 

3.4.1 The flow equation in the j —>■ oo Umit. 

We consider the flow of the Hamiltonian and an arbitrary observable O. Since H{i) determines 
U{£) one would expect a one-way coupling between the equations. It is assumed that O is a 
hermitian operator constructed in terms of Hq and Hi. Furthermore {z\0\ z) must be a rational 
function of j, which ensures that when taking derivatives of g{x,y,i) no additional factors of j 
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are generated. First we summarize the equations concerned 

H{i) = f + /^°AFo + f^AHi + i/ii {AHo, AH,} + ^f'AH^ + ^f^AHf + • • (3.13) 
0{i) = g + 5^°Ai7o + 5°'AFi + ^g'^ {AHo, AH,} + ^g^'^AH^ + ^g'^^AHf + ■■■ (3.14) 



dHji) 
dt 

doje) 

dl 



[[Ho, Hii)] , H{1)] , H{0) =Ho + \Hi, (3.15) 

[[Ho,H{e)],o{e)], 0(0) = 0, (3.16) 



where f{{Ho) , (Hi) , 0) = (Hq) + A (Hi) and g{{Ho) , (Hi) , 0) is just {z \0\ z) up to leading order 
in J. Next we substitute the expansion of H{t) into the flow equation (|3.15|) and take the 
expectation value on both sides with respect to the coherent state. Arguing as before we identify 
the leading order term on the left as being df /di. A general term on the right is of the form 

jinjkl ^[[^^^^ p^^^( . ^ p^^^(^ ^ ^3^^7) 

By transforming to scaleless variables and using the result of AppendixOit is seen to be at most 
of order 0(j2~*+L(i+*)/2j where t = i + n + k + l. Since the t = and t = 1 terms are zero 
(they involve commutators of scalars) the leading order contributions come from the t = 2, 3 
terms. These are exactly the terms found by considering the expansion of H{t) up to second 
order in the fluctuations. The expectation values of the double commutators may be calculated 
and expressed as functions of x and y using the method outlined in AppendixEl It is found that, 
due to cancellations, only four of the potential flfteen terms make leading order contributions. 
The corresponding expectation values are, to leading order 

([[Ai/o,Ai/i],Aifo]) = li{x,y)j (3.18) 

{[[AHo,AHi],AHi]) = 72(x,y)j (3.19) 

{[[AHo,{AHo,AHi}],AHi]) = ri{x,y)f (3.20) 

{[[AHo,AHf],AHo]) = -j3{x,y)f (3.21) 

where 71 = —4?/, 72 = 2x(l — x'^) and 73 = —2 + 4a;^ — 2x^ + 8y^. When keeping only the leading 
order terms on both sides the flow equation becomes 

% = lif'^'r + 12}'' f'' + l^sP'f' - llsf'f'', (3.22) 
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where superscripts denote derivatives to the rescaled averages x and y. For further discussion it 
is convenient to use the variables {x, q = cos(26')) G [—1, 1] x [—1, 1], where q and y are related 
hy q = 2y/(l - x^). The square domain of these variables also simplifies the numerical solution 
significantly. We note that the arguments above can be applied, completely unchanged, to the 
flow equation of 0{£) as well. We obtain, now for both H{i) and 0{i), the coupled set 

di \ dq dqdx dx dq^ j dq dx 

dg ^\ f ^9 dgd'^f\ ^ df dg (3 24) 

di \ dq dqdx dx dq^ ] dq dx 



Note that, in contrast to the equation for /, the equation for g \s a. linear equation that can be 
solved once / has been obtained from (|3.23j) . 

In Section 13.31 it was mentioned that H{() retains its band diagonal structure during flow, 
which means that Hi only appears linearly in the representation of H{£). This implies that 
/{{Hq) , {Hi) ,i) should be linear in {Hi), or, in the new variables, linear in q. When the form 
f{x,q,£) = nQ{x,i) + ni{x,i) q is substituted into 1)3. 23(1 this is indeed seen to be the case, and 
we obtain a remarkably simple set of coupled PDE's for no and ni: 

driQ dni 
-dl = -'"^5^' 

dni dno , . 

-di = -'"^ ^^-^'^ 

The initial conditions are 

Aw A(l — x'^) 

no{x,0) = x and ni{x,0) = — = '-. (3.26) 

q 2 

3.4.2 Finding the spectrum 

From / = jf it is possible to obtain the entire spectrum of H. For this discussion it is 
convenient to use the original {{Hq) , {Hi)) variables. Recall that in the i ^ oo limit H{i) flows 
toward a diagonal form and that the eigenvalues appear on the diagonal in the same order as in 
Hq, i.e. increasing from top to bottom. This implies that the n*^ eigenvalue of H is given by 

En = {-j + n\Hie = ^)\-j + n), (3.27) 



where Eq corresponds to the ground state energy. As H{i) flows towards a diagonal form the 
terms of expansion ()3.13|) containing Hi will disappear and eventually H{oo) and / will become 
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functions only of Hq and (Hq) respectively. The eigenvalues of H are given by 

En = f{{Ho) = -j + n,0,i = oo), (3.28) 

where (-f^i) has been set to zero. This can be understood in two ways. Looking at equations 
1)2. 4|) and ()2.1|) we see that the functional dependence of H{i) on Hq is the same as that of f 
(now called /) on (^^o)- Since taking the expectation value of H{oo) with respect to \—j + n) 
is equivalent to substituting —j + n for Hq, the result follows. Alternatively, consider equation 
(|2.11j) and note that setting {Hq) = —j + n makes the {n + 1)*^* diagonal element of AHq zero. 
Since Hi does not appear we see that when taking the inner product with | —j + n) only the 
scalar term, /, will survive. 

3.4.3 Calculating expectation values 

Next we return to the arbitrary operator O introduced earlier. The aim is to calculate 
{En\0\En) where \En) is the (unknown) eigenstate of H = Hq + XHi corresponding to the 
energy En- We can formulate this calculation in terms of a flow equation by noting that 

{En\0\En) = {En\Uii)uHi)OUii)U^{£)\En) = {E^, i \0{£)\ E^, £) (3.29) 

where \En,£) = U'^{£) \En), 0{£) = W {£)0U{£) and U{£) is the unitary operator associated with 
the flow of H{£). This equation holds for all £, and particularly in the £ ^ oo limit. Since 
\En-, oo) = \ + n) it follows that 

{En \0\ En) = {-j + n |0(oo)| -j + n), (3.30) 

and we conclude that the expectation value of O in the \En) eigenstate of H is simply the 
(n + l)*'^ diagonal element of 0(oo) in the Hq basis. Furthermore {z \ 0{£) \ z) has the character 
of a generating function in the sense that 



dz*'^dz^ 



oc (-j + n|0(£)| -j + m), (3.31) 

Z = 2*=0 



thus knowing which is just {z |0(£)| z) up to leading order in j, is sufficient to obtain all the 
matrix elements of 0{t) with high accuracy. However, unless an analytic solution for g is known, 
we are limited to numerical calculations for the low lying states. In particular, the ground state 
expectation value is found by setting z = 0. 
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3.5 Moyal bracket method 

In Section [2. 21 the Moyal Bracket method was used to obtain a simple, but general, realization 
of the flow equation as a partial differential equation in the variables a, [5 and t. For the present 
case this equation becomes 

'B^^M. = p)^H{a, /5, l)]^ , H{a, /3, , (3.32) 

where [f/, V]^^ = iO {Up Va — Vp Ua) is the Moyal bracket to leading order in 9. Here 6 = 2t: /{N + 
1) = 27r/(2j + 1) and so an expansion in is indeed controlled by the size of the system. As was 
mentioned earlier the model specific information enters through the initial condition, which we 
need to construct in terms of a and (3. We will not perform this construction for the Hamiltonian 
directly, but rather follow the more general route of constructing an irreducible representation 
of su(2), in terms of which both the Hamiltonian and observables can be readily obtained. 

We wish to find three functions Jz{oi,(3), J+{a,(3) and J^(a,l3) which satisfy the su(2) com- 
mutation relations with respect to the Moyal bracket [•, -J^. Note that the Moyal formalism has 
allowed an essentially algebraic problem, that of constructing a specific representation, to be 
reduced to that of solving a set of differential equations. We also remind ourselves that, since we 
use the Moyal bracket in a first order approximation, our representation can only be expected 
to be correct up to this same order. 

We begin the construction by making the following ansatz 

J+ = e'^f(a), J_ = e-'^f{a) and = p{a) (3.33) 

which clearly requires the representation to be unitary. This ansatz is based on the interpre- 
tation oi h = e^^ as a ladder operator which connects states labelled by the eigenvalues of 
g = e*". Substituting these forms into the required commutation relations [Jz, J±]^, = iJ± and 
[J-f, J_]^ = 2Jz produces the set 

e^p{a) = 1 and - e—f{a) = 2p{a) (3.34) 
da da 

which is easily solved to obtain 

2 

a o a a 

p{a) = -^ + ai and f (a) = --^ - 2ai- + a2. (3.35) 



Here ai and 02 are integration constants that we fix by requiring that the second order Casimir 
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operator assumes a constant value corresponding to the j-irrep: 

= j{j + 1) = ^{J+*J- + J-* J+) + Jz*Jz = a2 + al (3.36) 

We can satisfy this constraint to leading order in j by setting 02 = and ai = —j. Finally we 
arrive at 

Before continuing we define some more suitable variables. Consider the operator g = —i\og{g)/9 
which, according to equation ()2.17|) . has eigenvalues 0, 1, ... , 2j and is represented by a/9. This 
suggests that the natural domain of a/ 9 is [0, 2j]. We will use the scaleless variable x = a/ {j9) — 
1 G [—1, 1] in what follows. The representation now becomes 

J+ = je^f^Vl-x^, J- = je~'^Vl-x^ and Jz = jx. (3.38) 

We remark that this construction is by no means unique, although it is expected that other 
examples are related to this one by a similarity transformation. For example, complex non- 
unitary representation such as 

J+ = ie''^{jix-l)-l), J_ = ie~*^(j(a; + 1) + 1) and = jx. (3.39) 

may be constructed. In fact, this is an exact representation to all orders in 9, since x (and thus 
a) only occurs linearly and the Moyal bracket truncates after the first derivative. This illustrates 
the important role of the association between the eigenstates of g and the particular basis of the 
Hilbert space. For our purposes the unitary representation 1)3. 38(1 will be sufficient, although it 
is straightforward to check that the same results can be obtained using (|3.39j) . 

By substituting the representation (|3.38j) into the Lipkin model Hamiltonian, and being careful 
to use the Moyal product in calculating the squares of J±, we obtain the initial condition, to 
leading order, as: 

H{x,f3,0) = ja; + y(l-x2)cos(2/3). (3.40) 

The flow equation reads 

dH 1 

-Q^ = -T [HppHr^ — HpHpx) ■ (3-41) 

Again it is expected that the band diagonality of H(i) will be manifested as a constraint on the 
form of the solutions of the flow equation. We note that scattering between states of which the 
spin projection differ by two is associated with the cos(2/?) term also appearing in the initial 
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condition. Motivated by this we try the form H{x,l3,i) = j{no{x,i) + ni(x, £) cos(2/3)). Note 
that a factor of j, responsible for the extensivity of the Hamiltonian, has been factored out. 
Upon substituting this form into the flow equation we obtain 

dm . dm dn, dno 

which agrees with the equations of (|3.25j) in the previous section. For this form the flow of an 
observable is given by 

do , ^.dnidO , , 50 , , , 

_ = 2,,„(23)___4c„s(2;3)„,-. (3.43) 

We have seen that the governing equations obtained using the Moyal bracket method agrees 
with those of the previous section, although the interpretation of the constituents do differ. 
Next we turn to the matter of extracting the spectrum and expectation values from the so- 
lutions of these equations. As this procedure is largely similar to that of the previous sec- 
tion we will remain brief. First, note that through the representation (|3.39j) we may consider 
no(x,oo), and thus H{oo) itself, to be a function only of Jz- We arrive at the familiar result 
that En = jnQ{x = — 1 + n/j) for n = 0, 1, . . . , 2j. Expectation values are given by the diagonal 
elements of the flowed hermitian observable 0{x, l3,oo). The forms of J±{a,(3) suggest that in 
general 0{£) may be written as 

oo 

0{x, P,e) = Y^ fn (x, £) cos(2n/3). (3.44) 

n=0 

where cos(2n/3) corresponds to an off-diagonal term proportional to J^" -|- J^". (We only expect 
even powers since there is no mixing between the odd and even subspaces.) Importantly the 
/o(3;,oo) term contains the desired information about the diagonal entries. We can isolate /o by 
integrating over /3, which projects out the off-diagonal terms. Thus, in summary, 

{En\0\En) = {-j + n\0{(x>)\-j + n) 
= fo{x = -I + n/j, e = oo) 

= —[ dl3 0{x = -l + n/j,P,i = oo). (3.45) 
2vr Jo 

3.5.1 Probing the structure of eigenstates 

The ability to calculate expectation values for a large class of operators enables us to probe 
the structure of the eigenstates in a variety of ways. For this purpose we consider diagonal op- 
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erators of the form 0{f) = ^m=-j fi''^/j) 1"^) (^-1 where / is a smooth, j-independent function 
defined on [—1,1]. Since 0{f) = f{Jz/j), the initial condition for the flow equation is simply 
0{f,a,(3,e = 0) = /(a). Suppose = • Cm \m) is the eigenstate under consideration, 

in which case 

j j 

{E^\0{f)\E^)= J2 f{m/j)\cm\^= J2 fim/j)\{m\En)\^. (3.46) 

m=~j m=—j 

This leads to the interpretation of the expectation value as the average of the expansion coeffi- 
cients squared |(m|ii^rt)|^, weighted by the function f{x). Ideally we would like to choose 0{f) 
as the projection operator onto some basis state |m), as this would provide the most direct 
way of calculating the contributions of individual states. However, the projection operator does 
not fall within the class of operators corresponding to smooth initial conditions since there is 
no continuous function /(x) for which 0{f) = \m) {m\ in the large j limit. Instead we choose 
f{x,x) = exp(— 7(x — x)^) where x G [—1, 1] and 7 >> 1. This weight function focuses on the 
contribution of those basis states \m) for which m/ j lies in a narrow region centered around x. 
Considered as a function of x, it is expected that (0(/(x,x))) would approximate up to 

a constant factor, provided that the latter varies slowly on the scale of l/-^/7 in the region of 
X. Clearly 7 controls the accuracy of this method, and also determines the scale on which the 
structure of the eigenstates can be resolved. 

3.6 Numerical results 

In this section we analyze the results obtained by solving the equations derived in the previous 
two sections. We will mainly use the terminology of the Moyal bracket approach in what follows. 
See |23j for a treatment in terms of the expansion method. For compactness the £ argument of 
functions will occasionally be suppressed, in which case i = co should be assumed. 

3.6.1 Spectrum and expectation values 

First, we consider some structural properties of no(^) and ni{tj. It is known that the ma- 
trix elements of the Lipkin Hamiltonian possess the symmetry {m\H\m) = — {—m\H\ —m), 
which implies that the spectrum is anti-symmetric around Ej = 0. This symmetry is respected 
by the flow equation and manifests itself through an invariance of equations 1)3.42(1 under the 
substitutions uq — > —no, ni rii and x —x. This implies that no(x,£) = — no(— x,£) and 
ni{x,i) = ni{—x,i), and so we may restrict ourselves to the interval [—1,0] in the x dimension. 
This will be seen to correspond to the negative half of the spectrum, from which the entire 



3. The Lipkin Model 



29 




-1 -O.B -0.6 -O.-l -Q.2 ^ ZOO 4D0 500 

(c) (d) 

Figure 3.2: (a) and (b) The solutions to equation ()3.42p at ^ 100. The dots in (a) 
represent the exact spectrum for j = 1000. For clarity only every twelfth eigenvalue 
is shown, (c) The derivative of no with respect to the state label x. This is related 
to the gaps in the spectrum by equation ()3.47p . (d) The relative error in the first 
five states as a function of j for A = 0.5 

spectrum can easily be obtained. We also note that ni(ibl,0) = 0, and that this remains the 
case at finite i since dni/di is proportional to ni. Applying the same argument to no seems to 
suggest that no flow occurs for no at x = —1 and that no(— 1,^) = —1 for all i. This conclu- 
sion is, however, incorrect since it is found that in the second phase dni/dx develops a square 
root singularity at x = —1, which allows no(— 1,^) to flow away from —1. Numerically this can 
be handled easily by solving for m{x,i) = n1{x,i), instead of ni. We made use of the well 
established fourth order Runge-Kutta method jHl] to integrate the PDE's to sufficiently large 
values. 

It was shown that the eigenvalues of H are given by En = jno{—l + n/j). Furthermore, for 
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Figure 3.3: Energy levels and gaps as functions of the coupling constant, (a) The 
ground state energy, (b) The gap between the first excited state and ground state, 
(c) The gap between the second excited state and ground state, (d) The gap between 
^501 and S499 for j = 1000. 



sufficiently large j it holds that 
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dnc] 



dx 



(3.47) 



-l+n/j 



Figures 13.21 (a) and (b) show no and log(ni) at ^ ~ 100 for three values of A. At this point in 
the flow ni, which represents the off-diagonal part of H(£), is already of the order of lO"^'^, and 
may be neglected completely. As a comparison with exact results we calculated the spectrum 
for j = 1000 using direct diagonalization and plotted the pairs (x = —1 + n / j , E'^^'^^ / j) for 
n = 0, . . . , J as dots. We observe an excellent correspondence for all states and in both phases, 
with an average error of about 0.05%. This small discrepancy can be attributed to numerical 
errors and finite size effects, since we are comparing exact results obtained at finite j with those 
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Figure 3.4: (a) The order parameter Q = 1 + ( J^) /j as a function of A. (b) (J^) If' in 
the ground state as a function of A. (c) Expectation values of as functions of the 
state label x for different values of the coupling, (d) (J^) + 1 and ((AJ^)^) in I-E250) 
for 2 = 1000 as a function of A. The sharp local minimum occur at the point where 
x(A) = -0.5. 



of the flow equation which was derived in the j — > 00 limit. This is again illustrated in Figure 
13.21 (d) which shows the relative error in the first five eigenvalues as a functions of j for A = 0.5. 
The lines fall almost exactly on one another, so no legend is given. The log-log inset shows a set 
of straight lines with gradients equal to one, clearly illustrating the 1/j behaviour of the errors. 

Next we consider the dependence of the eigenvalues on the coupling A. Figure 13.31 (a) shows 
the ground state energy as a function of A together with the exact result for j = 1000. We see 
that Eq is fixed at — j in the first phase and begins to decrease linearly with A at large coupling. 
The phase transition in the Lipkin model is known to be of second order, and is characterized by 
a large number of avoided level crossings (SHI- This brings about complex non-analytic behaviour 
in the gaps between energy levels. The non-perturbative treatment of the flow equation enables 
us to reproduce much of this behaviour correctly. As examples we display the gaps E\ — Eq, 



3. The Lipkin Model 



32 




(a) 





Second 


Phase 


Pi rat Phaae 








S wond 


Phase 



(b) 



S 10 12 



-1 




(c) 



„ -! 

5 

= -s 




(d) 



Figure 3.5: (a) The negative half of the spectrum as a function of A. Every fourteenth 
eigenvalue is shown, (b) x(A) versus A. The relevant points in both the positive and 
negative halves of the spectrum are shown, (c) Energies of pairs of consecutive 
states of differing parity as function of A. Every third pair of the first 300 states are 
shown, (d) A log-log plot at A = 1 of no(x) + 1 versus x + 1, which is equivalent to En 
versus n. The slope of the linear fit is 1.33207; within 0.1% of 4/3. 



E2 — Eq and E^Qi — -E500 in Figure 13.31 One quantity that is not reproduced correctly are the 
gaps between successive states belonging to subspaces of differing parity. (Recall that subspaces 
of odd and even spin projection (parity) are not mixed by the Hamiltonian, as was shown in 
Section [3.21 ) As a definite case consider the gap Ei — Eq which is known to vanish like 1/j in the 
second phase, while the flow equation produces a gap which grows linearly with A. This can be 
attributed to the fact that the gaps are not extensive quantities and so they depend on higher 
order corrections in 1/j, which were neglected in our derivation. 

The function ni{x,i) found by solving equations (|3.42j) can now be substituted into equation 
(|3.43j) to obtain the flow equation for a general observable O. On a technical note, we found it 
advantageous to use expansion (|3.44|) to write the flow equation for 0{x, (3, i) as a coupled set of 
1 + 1 dimensional equations for the fnS, rather than treating it as a general 2 + 1 dimensional 
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PDE. First we consider the flow of Jz, which corresponds to the initial condition 0{x, /3, 0) = jx. 
Figure 13.41 (a) shows the results obtained for the order parameter together with exact results 
for j = 1000. We again find excellent agreement in both phases. Other observables can be 
considered by a simple modification of the initial conditions. For example O{x,P,0) = (jx)^ 
produces the second moment of Jz in the ground state, as shown in Figure |H3] (b) . The solution 
to equation (|3.43)) provides us with expectation values corresponding to excited states as well, 
which we find by evaluating /o, as defined in (|3.44|) . at the values x = —1 + n/j for n = 1, 2, . . .. 
Figure |S3](c) shows /o for O = at different values of the coupling strength. 

The non-perturbative application of the fiow equation method to the Lipkin model has pro- 
vided some interesting new insights into aspects of the phase transition. We will investigate 
some of these next. First we fix some notation. In the large j limit we may label states with the 
continuous label x S [—1, 1] through the association x ^ In this way a state is labelled 

according to its fractional position in the spectrum, for example x = —1 always corresponds to 
the ground state, while x = —1/2 denotes the state lying one quarter way up the spectrum. Now 
consider the behaviour of no{x) and Uq^x) depicted in Figures 13.21 (a) and (c) respectively. We 
note that for A > 1 there always exists a value of < 0, denoted by x{X), where the derivative 
of no{x) very nearly vanishes. Furthermore, x{X) must be a point of inflection of uq^x) since 
the derivative of rao(x) to x cannot change sign. If dno(x)/d3; were to become negative the flow 
equation would become unstable and cause ni{x) to grow exponentially, contradicting the results 
of Section [1.2.21 concerning the form of H{oo). This explains why we only observe a plateau at 
x(A), and no more drastic behaviour. At A = 1 this happens precisely at the ground state, i.e. 
x(l) = —1, while in the first phase no such point exists. In 3ir it was found that the energies of 
the low-lying states obey the scaling law 

En+j (X n^l^/N^'^ = N{x + 1)^/^ (3.48) 

at A = 1. This can be confirmed using the flow equations by considering the behaviour of no(x) 
close to X = —1. Figure 1531 (d) shows a log- log plot of nQ{x) + 1 versus x + 1 together with a 
linear fit which reproduces the power of 4/3 to within 0.1%. Earlier a direct link was established 
between the gaps separating successive eigenvalues and the derivative of no. This suggests that 
x(A) corresponds to a point in the spectrum with a very high density of states, brought about 
by a large number of avoided level crossings. Interestingly this point always occurs at the same 
absolute energy, namely Ei^ij^^(^X))j = ~ji although this corresponds to increasingly highly excited 
energies relative to the ground state. We believe that in the thermodynamic limit, contrary to 
Figure l3?2l (c) . the derivative of no, and the corresponding gap, should vanish completely at this 
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point. This has been confirmed in 36^ where it was shown that the gap is given by 

where = 2j. That our solution does not reflect this can be attributed to numerics, as this 
concerns a single point which is effectively "invisible" to the finite discretization used in the 
numerical method. The off-diagonal part of the Hamiltonian, represented by rii, also exhibits 
striking behaviour at x(A). In fact, upon returning to the flow equation 1)3. 42j) we make the 
interesting observation that at the point where the derivative of no{x) vanishes, ni{x) is not 
forced to flow to zero, but may in fact attain a non-trivial fixed point value. This is shown 
in Figure 1221 (b) which shows log(ni(x)) as a function of x. One clearly sees a sharp peak at 
the point where the derivative of nQ{x) nearly vanishes. The peak only occurs for A > 1 and 
moves to the right as A is increased. This is also consistent with the known result |2] that an off- 
diagonal element rriij of H{£) decays roughly as exp(— (E'j-i — Ej-if'i) at large £. This suggests 
a connection between quantum phase transitions, the corresponding disappearance of an energy 
scale (gap) |35j in the thermodynamic limit and the absence of decoupling in the Hamiltonian, 
also in the thermodynamic limit. 

The occurrence of this point in the second phase lends itself to the interesting interpretation 
of a "quantum phase transition" at higher energies. Indeed, apart from possessing some notable 
properties itself, it separates regions of the spectrum with markedly different characteristics. It is 
well known that states alternate between odd and even parity as one moves up in the spectrum. 
In the first phase these odd-even pairs are separated by a finite gap. In the second phase one 
finds a degeneracy between these successive odd and even states developing below the x(A) point. 
In particular this implies a degenerate ground state with broken symmetry in the second phase, 
although, as we have seen, this description may be applied to all states below x(A). Figure 
(c) shows a subset of eigenvalues as a function of the coupling, clearly illustrating this coalescing 
of pairs. Turning to Figure ITU fd). which shows the expectation value (^E^ij^^-^j \ Jz \ as 
a function of the state label x, we see a sharp local minimum occurring at a;(A) which separates 
the two phases. At the point x{\) itself the corresponding state \E(ij^^(^x^-^j^ is characterized by 
sharp localization around the |j, — j) basis state, also noted by [Hni- This is clearly illustrated 
by Figure Em (d) which shows a pronounced decrease in the spread of the state l-E'j/4) in the 
basis occurring at the point where x{\) = —0.5. We end this section with two "phase diagrams" 
which we hope will further clarify the discussion above. Figure l331 (a) shows the A dependency 
of a subset of the negative, even eigenvalues. For A < 1 the eigenvalues are confined between 
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Figure 3.6: The square root of the ground state expectation value of 0{f{x,x')) at 
A = 1.3 and for a range of x' values. The solid line is the exact modulus of the ground 
state wavefunction in the Jz basis for j = 4000. All results have been normalized to 
have the same maximum value. 



and —j. As A increases first the ground state (shown in bold) and then the excited states 
begin to cross the E = —j phase boundary until eventually, in the large coupling limit, only the 
\Ej = 0) state retains its first phase character. For finite j successive eigenvalues show avoided 
level crossings on the phase boundary. In the thermodynamic limit one would, however, expect 
that successive eigenvalues will coalesce as they cross the phase boundary; signaling a vanishing 
gap. A similar diagram, based in the label x rather than the energy itself, is shown in Figure 
13.51 fb). Keeping in mind the symmetries En = — i?2j-n and | En) = — (£'2j-n \ Jz\ E2j-n) 
this discussion can easily be adapted to apply to the positive half of the spectrum as well. 

These findings, published in 122], have stimulated further investigation of these phenomenon 
in the contexts of exceptional points 24^ and semi-classical approximations ^36^. These have 
shed light on the origins of the high density of states and localization occurring at x(A), as well 
as the mechanism responsible for the degeneracy between states of differing parity. 

3.6.2 Structure of the eigenstates 

In Section 13.5.11 we outlined a strategy whereby the structure of the eigenstates could be 
probed using the expectation values of a class of specially constructed operators. In this manner 
it is possible to approximate the modulus of the expansion coefficients of the state in the Jz 
basis. Here we present the numeric results of this procedure using the operators 0{f{x, x)) with 
f{x, x) = exp(— 7(x — x)^). We will consider the ground state at A = 1.3 and hope to find that 
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{Eq \f{x,x)\ Eq) k, |(m = xj\EQ)\^. Figure ESI shows that this is indeed the case. As expected 
larger values of 7 produce a more accurate reproduction of the form of the absolute wavefunction. 

3.7 Solution in the local approximation 

The numeric treatment of the previous section, although very effective as a computational 
tool, obscures some of the more subtle properties of the flow equation. In this regard an analytic 
approach provides some valuable insights, and we consider one such treatment next. 

At A = the eigenstates of the Lipkin Hamiltonian are simply the Jz basis states. Of course, 
this is no longer the case at non-zero coupling, although we expect that, in the first phase at least, 
the eigenstates will remain localized in the spin basis. The unitary transformation diagonalising 
H will reflect this by mainly mixing states of which the spin projections differ only slightly. This 
is manifested in the flow equations through the locality of the evolution of no,i(x,^). By this 
we mean that the flow at a point x = xi is very weakly affected by the values of no,i(a:,£) at 
points far from xi. We first consider the low lying states which can be investigated by solving 
no,i(x,^) within a neighbourhood oi x = —1. Since x = Jz/j (from representation (|3.38j) ) this is 
equivalent to restricting the calculation to states for which {{Jz) + j) 1 2 « 1; an approximation 
commonly used in this context. We begin by linearizing the initial condition about the point 
X = — 1, which gives 

no(x,0) = -l + (x + l) and 0) = A(l + x). (3.50) 

The form of no(x,0) is chosen to coincide with the parametrization of H{i) we introduce next. 
From equation (|3.42|) it is clear that the linearity of no,i is preserved during flow, which allows 
for the simple parametrization of H{i) as 

H{x, p, e)/j = -1 + ao{£){l + x) + ai{e){l + x) cos(2/3), (3.51) 

where ao(0) = 1 and ai(0) = A. The flow of the coefficients is given by 

dao dai 

— = -aiai and — = -aioo (3.52) 

where I has been rescaled by a factor of four. These equations leave ao(^)^ ~ oi(^)^ = 1 — A^ 
invariant, and since ai(oo) = we conclude that H{oo) = — j + \/l — A^ (j + Jz)- This correctly 
predicts the characteristic square root behaviour of the gap Ei — Eq = Vl — A^, as well as the 
low lying spectrum within a harmonic approximation. 
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The natural question arising here is whether this approach can be extended to allow for the 
treatment of the highly excited states as well. Indeed, if the flow at a point c G [—1,1] exhibits 
this locality property it seems reasonable that a local solution about c would provide a good 
result for the corresponding energy eigenvalue ^^(c+i). We will show how this local solution can 
be found analytically, and that by combining the solutions at various points an approximation 
for no{x,oo) may be constructed. The first step is again the linearization of H{0), now about 
an arbitrary point c G [—1, 1]. The linearity of H{£) in (x — c) is conserved during flow, allowing 
for the parametrization 

H{e)/j = ao{e) + ai{l) {x-c) + [a2{e) + a^il) {x - c)] cos(2/3), (3.53) 

where the flow of the coefficients are given by 

dao dai 2 da2 daa /o tr^\ 

^ = -di -dl = -M ^^-^^^ 

The initial conditions are ao(0) = c, ai(0) = 1, 02(0) = A(l — c^)/2 and 03(0) = — cA. These 



equations leave I = a\ — a\ invariant, and since 03(00) = we conclude that ai(oo) = Vl — c^A^. 
The presence of c in this equation is significant, as it indicates that the validity of the local 
approximation is a function both of the coupling and the specific point under consideration. The 
stability condition of the local solution is 

-\<c<\, (3.55) 

which holds for all c G [—1,1] in the first phase. In the second phase the stable domain is 
restricted to the interval [— 1/A, 1/A]. This is reminiscent of the discussion in the previous section 
concerning the transition appearing at higher energies in the second phase. It was seen that there 
exists points iba;(A) which separate regions of the spectrum for which the states exhibit either 
first or second phase behaviour. In fact, all the states lying within [—1/A, 1/A] have a first phase 
character, as x(A) is always found to be greater than 1/A. This suggests a connection between 
the phase structure, the properties of eigenstates, and the locality, or lack thereof, exhibited by 
the flow equation. 
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Figure 3.7: The results of the local approximation (LA) for (a) A = 1 and (b) A = 2 
together with the exact results for j = 1000. The dashed line corresponds to the first 
order solution given in ()3.58j) . Higher order local solutions, obtained numerically, 
are also shown. 

Continuing with the derivation, we find that the flow equations can be solved exactly: 

ao{£) = as (coth (aif + 02) — coth (02)) + c, ai{£) = ai coth (ai^ + 02) 

a2{i) = a3 cosech(Qi^ + 02), a3(^) = ai cosech(ai^ + 02). (3.56) 

The Oi's are integration constants that can be fixed as follows: 



01 



(00) = - C2A2 



ai = \/l — c^A^ 



ai(0) 



cosh(a2) 



-1 



03(0) 

"3 02(0) - 1 



ai 03(0) 



2c 



02 = arccosh 



-1 



"3 



cA 

- l)Vl - c^A^ 
2^ ' 



(3.57) 



The value of H{oo) /j at x = c is expected to provide a good approximation for the exact solution 
no(c, 00). We find that 



H(oo] 



ao(oo) = c + 



1 - c2) ( 1 - VI - C2A2 



2c 



(3.58) 



Figure ITTI shows this analytic result together with the exact eigenvalues. While it reproduces 
the spectrum of the first phase well, the local solution fares progressively worse at stronger 
coupling. This does not signal a lack of locality in the flow, but simply reflects the low order 
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to which the flow equation was solved. The inclusion of higher order terms ~ (x — c)" in the 
solution leads to much improved results, as is clear from the figure. Note that the point of 
breakdown for the fourth order solution in (b) has moved closer to the actual phase boundary 
at x{X = 2) —0.75 which separates regions of first and second phase character. 

In conclusion, we have seen that the local approximation provides non-perturbative results 
for states which possess a first phase character. In particular, the linear case can be solved 
analytically and reproduces the entire spectrum of the normal phase to good accuracy; the first 
such result that we are aware of. 



CHAPTER 4 
The Dicke Model 



4.1 Introduction 

Since its introduction in 1954 in the study of collective phenomenon in quantum optics, the 
Dicke model has received considerable attention in a wide range of fields, including that of 
quantum chaos and quantum phase transitions. We refer the reader to [SHI ^'^^ extensive list 
of references to these and related studies. Here we present a flow equation treatment of the Dicke 
model. Other boson and spin-boson problems have been considered in this context [391 Hj . and 
a special case of the Dicke model was treated perturbatively in [lUI- Vidal and Dusuel ,41^ have 
calculated finite-size scaling exponents for various quantities in the context of the Dicke model 
using flow equations. Compared to the Lipkin model this case presents a greater challenge to the 
flow equation approach, mainly due to the presence of two independent degrees of freedom. In 
order to keep the equations manageable we implement a two step procedure, first bringing the 
Hamiltonian into a block diagonal form before diagonalizing it completely. Unfortunately the 
complexity of the resulting PDE still prevents us from performing the first step exactly, and so 
a perturbative approach will be followed. 

We begin with an overview of the model and its structure relevant to the flow equation 
treatment. After introducing the necessary variables and representations we digress briefly to first 
treat a special case analytically before proceeding with the derivation of the full flow equation. 
Using this equation we will construct an effective form of Hamiltonian in which certain degrees 
of freedom have been decoupled. Finally the different approaches to diagonalizing this effective 
form are introduced and compared using numerical results. 

4.2 The model 

The Dicke model describes the interaction of N two-level atoms with a number of bosonic 
fields via a dipole interaction. For simplicity we will take N to be even; the odd case requires 
only minor modifications. Following [.3H|, we consider one such bosonic mode with frequency 
ui and coupling strength A. The corresponding bosonic creation and annihilation operators are 

and b. Associated with each atom are the spin-1/2 operators |s+\ s^'', s^''! which obey the 
standard su(2) commutation relations. All atoms have equal level splitting coq. The Dicke 
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Figure 4.1: A schematic representation adapted from [38j of the basis states 
{\m, n) : m = —j, ■ ■ ■ ,j; n = 0, 1, 2, . . .} for j = 2. Filled and hollow circles correspond 
to even and odd parity states respectively. Solid lines connect states belonging to 
the same Q-sector while dashed lines run between states connected by the J_|_6^ + J_6 
term of the Hamiltonian. 



Hamiltonian reads 

The I/^/N factor ensures that the Hamiltonian remains extensive when the bosonic mode is 
macroscopically occupied, i.e. when {b^b'j ~ N. By introducing the collective spin operators 
Jz = J2iLi ^z'^ ™d J± = J2^=i obtain the simplified form 

H = cjo J. + i^ib^b + ^ ( J+ + J_) (b^ + b) . (4.2) 

A natural basis for the Hilbert space TC is given by the eigenstates of the total collective spin 
operator = + + J^Jj^ together with and the boson number operator h = b^b: 

n = span {\j,m,n) : j = 0, . . . , iV/2; m = -j,...,j; n = 0, 1, 2, . . .} . (4.3) 

Here \j, m, n) = j{j + 1) \ j, m, n), Jz \ j, m,n) = m \ j, m, n) and h |j, m,n) = n \ j, m, n). Since 
J^] = the Hamiltonian does not mix different j-sectors of the total spin representation. 
The relevant sector for the low-lying states is j = N/2, and we fix j at this value throughout, 
dropping the j label in the basis states. We will focus on the resonant case where loq = loi = 1. 
As before, other cases simply correspond to different initial conditions for the equations we will 
derive. For these parameter values the model undergoes a quantum phase transition at Ac = 1 /2 
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from the normal to the so-called super-radiant phase. The second phase is characterized by the 
macroscopic occupation of the bosonic mode in the ground state, i.e. (n) ~ j. 

Next we consider the structure of H relevant to our application of the flow equations. Central 
to this discussion is the operator Q = Jz + n + j . With each distinct eigenvalue q = 0,1,2, .. . of 
Q we associate the corresponding subspace of eigenstates 



which we refer to as a Q-sector. The dimension of Qq increases linearly from a minimum of one 
at (7 = to a maximum value of 2j + 1 at q = 2j, from where it remains constant. A schematic 
representation of the basis states and Q-sectors appears in Figure 14.11 Grouping the terms in 
the Hamiltonian according to equation (|l.lfl|) leads to 



which makes it clear that H possesses a band block diagonal structure with respect to Q, as 
defined in Section ^2121 Since H only changes Q by two it leaves the subspaces corresponding 
to odd or even Q invariant. This implies a parity symmetry, similar to that of the Lipkin model, 
which may be compactly expressed in terms of the operator 11 = exp [inQ] as [H, IT] = 0. This 
symmetry gets broken in the second phase where the ground states of the odd and even sectors 
become degenerate. 

Early studies [441 145j of the model's thermodynamic properties were performed in the rotating 
wave approximation (RWA) which amounts to dropping the T2 + T_2 term. We will consider 
this case in more detail later on, as it reappears as the first order contribution to the solution of 
the fiow equation. 

4.3 The flow equation approach 

We have seen that the Hamiltonian contains interactions responsible for scattering both be- 
tween different Q-sectors and within a fixed sector. A complete diagonalization of H would 
require a transformation which eliminates both these interactions. Although the fiow equation 
is certainly capable of generating such a transformation in a single application, it is conceptu- 
ally clearer and technically much simpler to perform this diagonalization in two separate steps. 
As suggested in the previous section we will first eliminate interactions which connect different 




(4.4) 




(4.5) 
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Q-sectors through the structure preserving flow generated by r]{i) = [Q,H{£)]. As before this is 
done in the thermodynamic or large j hmit. A physical interpretation of Q is that of a counting 
operator for the elementary excitations or energy quanta present in the system. The flow will 
produce an effective Hamiltonian H^ff = H{oo) which conserves the number of these excita- 
tions. We expect there to be new terms, both diagonal and off-diagonal, appearing in H^ff- 
Apart from the provision that they conserve Q there is no obvious constraint on these newly 
generated interactions, and so no band (block) diagonal structure need be present within the 
Q-sectors. 

Having brought H into a block diagonal form there are two ways in which to proceed. We 
may construct, for large but flnite j, the matrices corresponding to each Q-sector and then apply 
direct diagonalization. Although still numerically intensive this is a greatly reduced problem 
compared to diagonalizing the original fully interacting Hamiltonian, as the submatrices are at 
most of size {2j + 1) x [2j + 1). Alternatively we may remain at inflnite j and apply the flow 
equations to each sector separately using, for example, the generator r/(£) = [Jz,H{iy\. For a 
flxed value of Q the problem now becomes effectively one dimensional as we may eliminate either 
the bosonic or spin degree of freedom in favor of Q, which acts as a scalar parameter. Both these 
methods will be demonstrated later on. 

Finally, a remark on the ordering of eigenvalues. When applying the flow equations to each 
sector in the second step we expect complete diagonalization and that the eigenvalues will appear 
on the diagonal in increasing order [^. For the flrst step, which ends in a block diagonal form, 
the situation is no longer so clear, as the ordering proof in „20 is only valid for the case of 
complete diagonalization. In short, it is generally unknown in which sector a certain eigenstate 
of H will be found when we diagonalize H^fj. In the flrst phase we expect the same locality 
encountered in the Lipkin model f Section 13. 7() and so the low-lying states should belong to the 
sectors with Q <^ j. In numerical investigations for small j values we have observed this ordering 
in both phases. In particular, the ground state is mapped to the single basis state of the Q = 
subspace. As this ordering is a non-perturbative phenomenon we do not expect to observe it in 
our perturbative treatment. 

4.4 Flow equations in the j — oo limit 

In the subsections that follow we derive the flow equation for the flrst step of the diagonal- 
ization process using the Moyal bracket method. We begin by introducing the relevant variables 
and representations. 
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4.4.1 Variables and representations 

To account for the model's two independent degrees of freedom we introduce two pairs of 
operators g,h and g',h' as described in Section [2.21 These satisfy the exchange relations 

hg = e-'^gh and h' g' = e''^' g' h' (4.6) 

while operators coming from different pairs commute. The pair g, h is used to represent the 
spin degree of freedom through the representation constructed in Section 13.51 in terms of Xg = 
a/ije) - 1 G [-1,1] and /3: 

J+=je'^^l-xl J.=je-'^^l-xj and J,=jx,. (4.7) 

In similar fashion we wish to construct a representation for the boson algebra {6^, b} in terms of 
g' and h' . Two new issues, not encountered in the su(2) case, arise here. Firstly, it is well known 
that only infinite dimensional representations of the boson algebra exist, whereas g' and h' are 
finite dimensional. The second point concerns the scale we should associate with the bosonic 
operators. If b^ and b were naively assumed to be scaleless with respect to j (or N), only the 
Jz term in the Hamiltonian would survive when working to leading order in j. We will address 
these issues in the context of our proposed approximate representation 

b = h'^^, b'^ = ^h' and bh = n = g (4.8) 

where g' = —i log(g') /O' = diag(0, 1,2,... , — 1) and D' is the dimension of the space. Denoting 
the eigenstates of ^' = n by |n) n = 0, 1, . . . , D' — 1 we see that b\n) = y/n\n — 1) and 6^ \n) = 
\/n + \\n + 1) for < n < D' — 1. (See equations (|2.17|) and (|2.18|) .) The creation operator 6^ 
maps the highest state \D' — 1) to zero, since h' \D' — 1) = |0). This amounts to a truncation 
of the boson Fock-space, and the operators in (|4.8() agree with the truncated forms of the exact 
infinite-dimensional operators. Proceeding as before, we treat g' and h' as scalars and define a' 
and /?' through g' = e*" and h! = e*^ . The representation now becomes 

b = e-'^' ^/a'/O' + 1, ftt = e'f^'^a'/e' and bh = a'/e' (4.9) 

which satisfies, to leading order in 6' , the desired commutation relation with respect to the 
bosonic Moyal bracket ^^]^ = 1- For the moment we set D' = uj + 1, and return to the role 
of u later. The second issue concerns the scale of j to be associated with b and b^ . This is really 
a question of which values of a' /9' are relevant to the low energy physics of the system. In the 
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Dicke model it is known |SH] that ground state occupation of the bosonic mode is microscopic, 
i.e. (n) 0{j^), in the first phase and macroscopic, i.e. (n) in the second. We 

can describe both these cases by considering b^b = a' /9' oc j. A natural choice of variables is 
Xb = a'/{j9') G [0,u] which is scaleless and analogous to Xg of the su(2) case. The j-scale now 
becomes explicit in both the representation 

b= Vje-'^'y^xh + l/j, b^ = y^e'^'^b and bh = jxb (4.10) 

and bosonic Moyal bracket 

[U{xb, /?'), V{xb, 13')] ^ = i {Up,V,, - Vp,U,,) /j. (4.11) 

We also observe that the dimension D' = uj + 1, which controls the Fock-space cutoff, does not 
appear explicitly. For practical purposes we may safely assume u to be much larger than the 
interval of x^ under consideration. 

Combining the two representations we obtain the initial condition to leading order in j as 

H{xs,P,Xb,(3',i = 0) =j \xs + xb + A ^2x6(1 - x2) cos {(3 - /?') + A ^2x^(1 - x2) cos {(3 + /?')] , 

(4.12) 

where the terms are ordered as in equation (|4.5j) . 

4.4.2 Solution in the local approximation 

Before proceeding with the derivation of the general flow equation, which will involve signifi- 
cant behind-the-scenes numeric and symbolic computation, let us consider a case amenable to an 
analytic treatment. Using the same locality argument put forth for the Lipkin model in Section 
13. 71 we assert that for the low-lying states only the flow within a small region around = — 1 is 
relevant, and that the evolution of H in this region is governed by its local properties. This is 
equivalent to the assumption that ( -|- j) ~ for the low-lying states. A similar approximation 
is made in the bosonization treatment of |SH1) the results of which will be reproduced here using 



the flow equations. In practice we implement this approximation by simply replacing \J\ — x^ 
by \pl^\ + Xs in the initial condition (|4.12|) : the result of a Taylor expansion to leading order in 
Xs + 1. First we consider the flow generated by ??(£) = [Q, -ff(^)]. As with the Lipkin model the 
local approximation leads to a very simple form for H{i), with only two new terms, proportional 
to 6^ -)- [b^)"^ and + J^, being generated. As required the band block diagonal structure of H 
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is conserved. The flowing Hamiltonian may be parametrized as 



H{£)/j = -1 + aoW (xs + Xb + l) + 2 ai(£) ^^^^(l + x,) cos (/? - /?') 



(4.13) 



Xfo COS (2/?') + {xg + 1) cos (2/?)) 



where aj(£) i = 0, 1,2,3 are scalar coefficients. In operator language this amount to 



H( 



-j + ao{£) (j + Jz + n) + ai{£) 



J+h + J_6T 



(4.14) 



+a2{£) J+b^ + J-b /^/j + a^{£) {{b^ f + b'') + (JI + J1) / j /2. 



Observables linear in and h may be parametrized in a similar fashion. We will consider 
= Q — j = Jz + n for which 0{£) is of the form ()4.13|) with the flowing coefficients denoted by 
ai{£) for i = 0,1,2,3. 

Substituting these forms into the flow equations (|4.31|) and (|4.35|) and then matching the 
coefficients on both sides yield 

dao 2 2 dai ^ da2 da^ 

-^ = -a2-«3' -^ = -20203, = -00^2 - ai03, — = -0102 - aoas (4.15) 

and 



dao _ _ dai 

— = -a202 - 0303, — 
ai dl 



da2 _ _ da^ 
-0203-0302, -7^ = -0002 - 0103, -j-r = -0102 - 0003. 
d£ dl 



(4.16) 

The initial conditions for H{£) and 0{£) are oo(0) = 1, ai(0) = 02(0) = A, 03(0) = and oq = 1, 
oi = 02 = 03 = respectively. 

First we consider H{£). In the ^ — > 00 limit 02 and 03 is expected to vanish, while oq and 
oi will assume new, renormalized values. The initial and final values of the Oj's may be related 
using the fiow invariants /i = Oq + o^ — o| — 03 and I2 = oqOi — 0203. By combining these with 
the initial conditions we find that 



Oo(c!o) + Oi(oo) = 1 and 00(00)01(00) = A. 



(4.17) 



The renormalized values of oq and oi correspond to a point where the circle and hyperbole 
described by these equations intersect. For A > 0.5 no such point exists, and the solution breaks 
down. For A G [0,0.5] the fixed point is found to be 



Oq 



;oo) = (Vl + 2A + Vl - 2a) /2 and oi(oo) = (Vl + 2A - Vl - 2a) /2. (4.18) 
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(a) (b) 

Figure 4.2: The flow of {ao{i),ai{£)) for (a) A = 0.45 and (b) A = 0.6. The dashed Hnes 
represent the circle and hyperbole described by equations (|4.17|) . In (b) the flow is 
only shown up to some finite value of i, as no convergence takes place. 



This is illustrated in Figure which shows the flow of {ao{i), ai{i)) for A = 0.45 and A = 0.6. 
Next we consider the flow of 0{i). To fix the final values of the coefficients we require four 
invariants, some of which must involve the coefficients of H(£). One such set is 

7i = Og + «! — — a|, I2 = aooi — 0203, I3 = aooo + aioi - 0202 — 0303 

and = ooai + aioo — 0203 — a3a2- (4.19) 

Combining 73 and I4 with the initial conditions lead to 

oo(oo)ao(oo) + ai(cxD)ai(cxD) = ao(0) = 1 

ai(cxD)ao(cxD) + ao(cxD)ai(cxD) = ai(0) = A (4.20) 



and so 



, , 1 / 1 + A 1- A 
00(00) = - , + 



,,1/1+A 1-A, 
ai(oo) = - -j== - -j== . (4.21) 



V1 + 2A VI - 2A 
1+A 1-A 
2 VVl + 2A VI - 2A 

The values of 02(00) and 0,3(00) can now be solved for using Ii and /2- Note that for the purposes 
of calculating expectation values with respect to the eigenstates of H only ao(oo) and ai(oo) are 
relevant. This follows from the observation that 02 and 03 correspond to terms that change 
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Q by two, whereas the eigenstates of Heff are also eigenstates of Q. For simpHcity we drop 
these "between-sector" scattering terms in what follows, keeping only those terms relevant to 
the calculation of expectation values. 

The successful construction of the effective Q-preserving Hamiltonian i?e// = H{oo) and the 
observable Oe// = 0(oo) marks the end of the first step. We now apply the flow equation a 
second time, using the generator r]{£) = [Jz, Heff{i)]. When restricted to a single Q-sector 
is clearly non-degenerate, and so we expect complete diagonalization in the final Hamiltonian. 
The local approximation again allows for a simple parametrization 



Heff{i)/j = -1 + bo{e){xs + 1) + bi{e)xb + 262^7^6(1 + COS (/? - /?'), (4.22) 

where bo{0) = 6i(0) = ao(oo) and 62 = ai(oo). The flow equations for the coefficients are 

^^^l ^^-^l a.d f ^-(^.-W^ (4.23, 

which leave = 60 + ^1 and l'^ = 69 + ^1 + ^^i invariant. Combining these with the initial 
conditions and using the fact that 62(00) = leads to 

60(00) = ao(oo) + ai(oo) = Vl + 2A 



61(00) = ao(oo) - ai(oo) = Vl - 2A. (4.24) 

Turning to Oeff{£) we find that its parametrized form is again similar to that of Heff{£), and 
we denote the flowing coefficients by bi{£). The flow equations are found to be 

^ = 26262, ^ = -26262 and ^ = -(60-61)62 (4.25) 
which leave J( = 60 + 61 and I2 = 6960 + 6161 + 26262 invariant. Proceeding as before we find 



that 60(00) = (1 + A)/\/l + 2A and 61(00) = (1 — A)/\/l — 2A. Since the eigenstates of i?e//(oo) 
are also eigenstates of Jz we drop the 62 term. 

In summary, we have seen that within the local approximation the two step diagonalization 
procedure can be performed exactly. The final Hamiltonian and transformed observable O = 
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Jz + n is 

H^ff (oo) = -j + {j + Jz) Vl + 2A + nVl - 2A 

o,„M = -; + 0- + J=)^ + ^^ + ... (4.20) 

which agrees with the result of j38j obtained using bosonization. This clearly constitutes a 
harmonic approximation of the spectrum in terms of two oscillators with frequencies \/l + 2X 
and \/l — 2A. Only the parts of Oeff which commute with both Q and Jz are shown. Note that 
the gap between the ground state \Eq) = \—j,0) and first excited state \Ei) = exhibits 
the characteristic square root behaviour A = \/l — 2A. 

4.4.3 Flow equation for the Hamiltonian 

We proceed with the derivation of the general flow equation for the Hamiltonian, beginning 
with a more in depth study of the structure of H{t). In Section fl. 2. 31 it was shown that the flow 
generated by r]{t) = [Q, will preserve the band block diagonal structure present in H{0), 
thus restricting the terms in H{£) to those changing Q by either zero or two^ 

oo 

H{e) = fo{Jz,n,i)+fi{Jz,n,£) (^J+b^ + J_6)+ fkAJz,n,<i) ((J+)'(6y + {J-)Hb^y) ■ 

k=l 

le{k,k±2},k+l>2 

(4.27) 

A simpler and much more convenient form can be found by introducing the operators Tq = 
J+b + J-b^ and 

^2 = {Ts = Jl + Jl, n = {h^ f + = J+b^ + J-b] . 

We interpret Tq as the fundamental Q-preserving interaction while the elements of J-2 constitute 
the fundamental interactions which change Q by two. As shown in Appendix^ we can rewrite 
all the operators appearing in 1)4. 27|) in terms op Jz, n, Tq and the elements of J-2 which only 
appear linearly. For example: 

J+63 + J_(6t)3 = Tgn-hTsb 
(J+)V + (J_)2(6t)2 = r2 + 2n(jf-j2) 

(J+)V + (J-)^(6t)^ = {{j^^-f)n + T^)T, + {j!-f)TqT,,. (4.28) 

^In this equation, and in some that foUow, we treat operators as commuting scalars. Since reordering can only 
bring about 0{l/j) corrections this is sufficient for our purposes, and it simplifies the notation considerably. 
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This results in the form 

H{£) = j [Fo( J„ n, Tg,£) + Fi{J,,n, T„ e)Ts/f + F2{J,,n, T„ £)n/j + J,, n, T^, £)r,fe///2 

(4.29) 

The factor of j responsible for the extensivity of H has been factored out explicitly, and so each 
Fi is a scaleless function. The initial condition is given by 



Fo = Xs + Xb + ^Xq, Fi=F2 = and ^3 = ^ (4.30) 



where Xq = Tqj f'l'^ . We proceed by inserting this form into the flow equation 

— = [[Q,i7],,i/], (4.31) 

where [•, -J^ is the full Moyal bracket 

[[/, V\, = i {U^V^^ - VpU^^ + UfS'V^, - Vf,,U^,) /j. (4.32) 

The resulting algebra is easily handled using a symbolic processor such as Mathematica. Matching 
the coefficients of Tg, Tf, and Tsb on both sides of 1)4. 31() produces a set of equations of the form 

^ = Pt'''^tlu<n^l''' ^ = 0'l>2,3 (4.33) 

where we sum over repeated indices. The superscripts of the i^j's denote derivatives to the 
scaleless variables Xs = Jz/j, Xb = h/j and Xq = Tq/p^'^. The v,v' indices run over the set 
{0, Xs, Xh, Xq} where corresponds to no derivative being taken. Governing the flow of Fi is the 
16 X 16 matrix of which the entries are simple polynomials of the rescaled variables. 

The non-zero entries of these, typically sparse matrices appear in Appendix El Due to the 
obvious complexity of these non-perturbative equations we have been unsuccessful in finding an 
exact numerical solution. A perturbative solution can be found, and we present this in Section 
1131 

4.4.4 Flow equation for an observable 

Next we derive the flow equation for an observable O. This case is complicated by the lack 
of structure present in 0{l) at non-zero i. For example, we do not generally expect band block 
diagonality to be present, and so 0{i) may contain interactions connecting distant Q-sectors. In 
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general we are confronted by the form 



oo A rr 



A=0 p=0 
A even 



0W= E E^(A,P)(«^^'-'?;.^)7(^ (4.34) 



where T,^a,p) = {J+y{b^)^~P + (J_)p(6)^~p and T(o,o) = 1- The A label indicates the amount 
by which the term changes Q, and similar for p with respect to Jz- We define the index set 
Z = {{A,p) : A = 0, 1, 2, . . . ; p = 0, . . . , A}. Inserting this form into the flow equation 

^ = [[Q,Hl,Ol (4.35) 



yields the coupled set 



We sum over the repeated indices / G {0, 1, 2, 3}, /' G X and v, v' G {0, Xs,Xb, Xq\. is a 

16 X oo dimensional matrix containing polynomial functions of the scaleless variables. Note that 
at £ = CO only the i^(o,o) function is of interest for the calculation of expectation values. This 
follows from the observation that eigenvalues of H{oo) are also eigenvalues of Q, and that -^(o,o) 
is the only term leaving Q invariant. 

Again it is clear that a direct numerical approach is intractable, and we instead try to find 
perturbative solutions. 

4.5 Perturbative solutions 

We wish to construct solutions to equations (|4.33|) and (|4.36() of the forms 

Pmax Pmax 

H{£) = E X^H^^Ki) and 0{e) = E X^O'^p\£). (4.37) 

p=0 p=0 

The forms given in 1)4.29(1 and (|4.34|) are still valid, as they apply separately to each H^^^ and O^^^ . 
When working to finite order in A the functional dependence of H{i) and 0{i) on {x^, Xs,Xq} are 
constrained to be polynomial and of finite order. This also limits the type of interactions that 
can be generated in 0{£) since terms for which A = 2n is at least of order A", and so only F(^A,k) 
with A < [pmax/2j are relevant. This is true at £ = 0, provided that O ^ and [0,Q] = 0, 
and continues to hold at £ > since rj{i) is always at least of order A and contains only A = 2 
terms. As all the relevant functions are simple polynomials we may proceed by constructing a 
set of coupled differential equations for the scalar coefficients appearing in these polynomials as 
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functions of i. The extensive algebraic manipulations involved in this step can be automated for 
arbitrary Pmax using Mathematica. The resulting set of ordinary differential equations is then 
solved using the standard Runge-Kutta algorithm Solutions were obtained for H up to 

Pmax = 20 and for the observables Jz and n up to Pmax = 18. These transformed operators will 
be denoted by Hefj = H{oo), h^ff = n(oo) and {Jz)eff = '^2(00). We only state the results up 
to fourth order here: 



where, to aid interpretation, we temporarily abuse notation by writing Jz,n,Tq for the scaleless 
variables Xs,Xh,Xq. The rationality of the coefficients, apart from the \/2 factors, have been 
verified to very high numerical accuracy. In similar fashion the observables are given by 



jW(00)/j = Jz + ... 

ji'H^)/j = + ... 

4^\^)/j = {l-2nJz-Jl)/A + ... 

jf\oo)/j = (3n + 16J^)T,/(32V2) + ... 

)(cx))/i = (-18n - SlTg^ - 228 J;, - l2n^Jz + 282nj2 + 228^) /384 + . . . (4.39) 



fv- '{00)/ j = h + . . . 

n^^\oo)/j = 0+... 

n(2)(oo)/j = {i-2nJz-Jl)/A + ... 

n(3)(cx))/j = (5n + 16J^)rg/(32V2) + . . . 

h^^^ (00) /j = (I8n - 69r2 - 132 - 12n + 438n J^^ + 132 jf) /384 + . . . (4.40) 



where only the Q-preserving terms relevant to the calculation of expectation values are shown. 
Looking at the results for H(,ff we see the "in-sector" scattering interaction Tq appearing in the 



FW(oo)/j 

//«(CX))/J 

i7(2)(oo)/j 



n+ Jz 



Tq/V2 

( j2 + 2nJz - 1) /4 

-Tq {h + AJz) /{8V2) 

{7T^ + 20 -2h + 2,SnJl + 20Jf ) /64, 



(4.38) 



and 
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higher order terms. As the power to which Tq can occur is hmited by the perturbation order a 
band diagonal structure is present in each Q-sector with respect to Jz (or equivalently h). We 
expect Heff = H{oo) to provide a good description of the Dicke model at weak coupling, and 
possibly throughout the first phase. At strong coupling there is no guarantee of the accuracy 
or stability of H^jf, even when working to arbitrarily high orders in the perturbation. In fact, 
from the terms given in (|4.38|) two things become apparent which signal a breakdown of this 
approximation at A > 0.5. Firstly, the matrix element of the one-dimensional Q = subspace 
is found to be {—j,0\Heff \ — j,0) = —j for all A and orders of the perturbation considered. As 
is known from non-perturbative numerical studies the flow equation maps the ground state to 
precisely this \—j,0) state. However, it is also known that the ground state energy is only equal 
to — j in the first phase, and that it decreases linearly with A in the second [HHl) contradicting the 
prediction of H^jf. Secondly H^jf is found to become unstable at large coupling for perturbation 
orders higher than one. For example, consider the expectation value of H^ff to second order with 
respect to a simple variational state: {—j, n \Hf,ff \ — j, n) = —j + (l — A^/2) n. We note that for 
A > \/2 this expectation value is unbounded from below. Together with the variational principle 
this implies that the spectrum of -ffe// becomes unbounded from below at some A < \/2. For 



the fourth order case this instability occurs at a A < y 2(-v/TT — l)/5 ~ 0.963. We have observed 
that this point of breakdown continues to move closer to the critical point Ac = 0.5 as higher 
order corrections are included. This agrees with the general notion that a perturbation series, 
even when summed up completely, may produce divergent results beyond the series' radius of 
convergence. 

Finally we point out that to first order the effective Hamiltonian is 



which is equivalent to the RWA approximation in which the model was originally studied [441 145j . 
To this order we also observe that the Q-preserving parts of Jz and n are left unchanged by the 
transformation generated by the flow equation. 

4.6 Diagonalizing H^jf 

The flow equation treatment of the preceding sections has brought the Dicke Hamiltonian 
into a block diagonal form by eliminating interactions which connect different Q-sectors. What 
remains is to diagonalize Heff within each Q-sector at a time. 

One approach, valid for large but finite j, is to apply direct diagonalization to each Q-block 
(submatrix) of i^g/ / • This is done by first constructing the matrix representations of , n and Tq 





(4.41) 
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for the Q-sector under consideration. Finding the desired submatrix of H^jf is simply a matter 
of replacing Xg, Xb and Xq by the properly scaled matrix representations of Jz/j, n/j and Tq/j^^"^ 
respectively. One need not be concerned about the precise ordering of operators in products, as 
different orderings only bring about 0{l/j) corrections. The result of these substitutions is a 
matrix M, the size of which ranges from 1 x 1 to (2j + 1) x (2j + 1) depending on the value of 
Q. Generally M will not be Hermitian, and so we consider the symmetrized form (M + M'^)/2. 
Results of the subsequent numeric diagonalization appear in the next section. 

Alternatively we may remain in the large j limit and apply the flow equation again, this 
time non-perturbatively, to diagonalize each submatrix. Within each Q-sector the problem is 
effectively one-dimensional as we may eliminate either the bosonic or spin degree of freedom in 
favor of the scalar parameter Q. We choose to work with the spin degree of freedom, which has 
the natural basis 

SQ = {\m,Q) : m = -j,...,mm{j,Q - j)} . (4.42) 
By writing n = Q — — j and 

Tq = J+b + J.b^ = yjQ-Jz- j J- + - - j + 1 J+ 
the effective Hamiltonian can be rewritten purely in terms of spin operators. 

Let us clarify this last step. Since m and Q together fix the number of bosons there is no 
n label in Bq. In fact, we completely eliminate the bosonic degree of freedom through the 
correspondence 

J+h ^ y^Q - Jz - j + 1J+ and J^b^ ^ \/Q- Jz- jJ~ (4.43) 

where the operators on the right act on Bq and are equivalent to the operators on the left acting 
within a fixed Q-sector of the original tensor product space \m) ®\n). 

For concreteness we consider the fourth order case for the rest of the section. The results 
can be easily extended to higher orders. The form of the spin-only Hamiltonian is -f^e// = 
No{Jz) + Ni{Jz) (J+ + J-) + N2{Jz) {j\ + which clearly possesses a double band diagonal 
structure with respect to Jz- To preserve this form during flow we use the generator of Section 

r?(^) = Ni{Jz4) {J+ - J-) + N2{Jz,i) (J+ - J-) . (4.44) 
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For the numerical treatment of the flow equation it is convenient to transform to the variables 
of the su(2) representation in terms of which Heff{l) and r]{i) become 

Heffii) = j(no(x,£) + ni(x,£)cos(/3) + n2(x,£)cos(2/3)), (4.45) 
r]{i) = ij {ni{x,i)sm{(3) +n2ix,i)sm{2f3)) , (4.46) 

where we write x for Xg and the initial conditions are 

no(x,0) = g-l-(l-2(g-l)x + x2)AV4 + ((l + x)2(8x-3) + g(3-13j;2))AVl6, 



ni(x,0) = y/q-1- xVl - x^{8X - {q - I + 3x)X^)/\/32, (4.47) 
n2(x,0) = -7(g-l-x)(x2-l)A^/32. 

Here q = Q/j, which is now a continuous scalar parameter labelling the relevant sector. The 
flow of these functions are given by 

driQ dni dn2 

—— = —rii— 2n2— — , 

oi ox ox 

drii dni ( dnn dno \ , . . ^ 

-2n2^ -nA-^ + -^], (4.48) 



di dx \ dx dx 

— - -2n2 — 
d£ dx 

The corresponding flow equation for an observable 0{i) is 

do ( dni ^ dn2 A dO , ^ ^. dO , . . s 

' sm/3 + — ^sm2/3 — - (nicos/3 + n2Cos2/3) — , (4.49) 



dx dx J dj3 dx 

which we treat using a similar approach to that of equations (|3.44|) and (|3.45|) . In the £ ^ oo 
limit we expect ni and n2 to vanish, as they represent the off-diagonal parts of the flowing 
Hamiltonian. Using the correspondence between x and Jz we may consider no(cxD), and thus 
H(oo) itself, to be a function of Jz- It follows that the eigenvalues of the particular Q-sector are 
given by = f^o(2; = — 1 + n/j, i = oo) for n = 0, 1, . . . , min(2j, Q). 

There is a subtle issue regarding the domain of these functions that still needs to be addressed. 
Looking back at the definition of the Bq basis we see that the maximum value of the spin label m 
is min(j, Q—j) which is less than j when Q < 2j. Correspondingly we need to restrict the domain 
of the rii functions to reflect this. Since Jz oc x the relevant domain is x G [— l,min(l,g — 1)]. 
Indeed, outside this region the initial conditions are complex, causing the flow to become unstable 
as ^ ^ oo. 
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4.7 Numeric results 

During the course of this chapter we have encountered a number of different treatments of the 
Dicke model. In this section we present the results obtained using these methods and compare 
their ranges of applicability. First we summarize the different approaches: 

• Direct Diagonalization: By introducing a cut-off in the maximum value of Q we can 
numerically diagonalize the full Dicke Hamiltonian on a truncated Hilbert space. This 
generally leads to very large matrices, especially in the second phase where an extensive 
number of Q-sectors must be included in order to describe the low-lying states. For suf- 
ficiently large values of the cut-off this method does provide very accurate results for the 
low-lying states, and we will use these as a benchmark for the flow equation results. 

• Local Approximation: This analytic method was derived in Section r4.4.2l and constitutes 
a harmonic approximation of the first phase, low-lying spectrum in terms of two oscillators. 

• Direct Diagonalization of Heff- In Section 14.51 we derived the effective Q-preserving 
Hamiltonian H(.ff in a perturbative approximation. For finite j the submatrices comprising 
the block-diagonal structure of H(,ff may be constructed and diagonalized numerically. 
Since these submatrices are maximally of size (2j + 1) x (2j + 1) this is a greatly reduced 
problem compared to diagonalizing the original Hamiltonian in some large truncated space. 

• Flow Equation Diagonalization of H^ff. As described in Section 14.61 the individual 
Q-sectors may be diagonalized in the j ^ oo limit using a flow equation equipped with an 
appropriate generator. The continuous variable q = Q/j G [0,c>o) appearing in the initial 
condition labels the relevant Q-sector. 

• The Rotating Wave Approximation (RWA): Within this approximation the Dicke 
Hamiltonian is equivalent to the flrst order form of -f^e//- In this sense the RWA is already 
included in the two previous cases. However, due to its prevalence in the literature we will 
treat it separately reference case. 

Throughout this section only the subspace corresponding to even values of Q, as discussed in 
Section [4.21 will be considered. The first results we present are those obtained by diagonalizing 
Heff at finite j. The fiow equation treatment of the second step will be dealt with later on. 
Figure (a) shows the energy of the third excited state IS3) as a function of X at j = 100. Also 
shown are the exact values obtained using direct diagonalization, and the RWA result. Although 
the result of -f^e// is a marked improvement over that of the RWA, we still see 0{l/j) errors 
present at larger values of A. This simply the reflects the order to which the Moyal bracket was 
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(a) (b) 

Figure 4.3: The (a) third excited state and (b) the gap A£^3 as functions of A. i^e// 
was constructed for j = 100 up to perturbation order Pmax = 15. 



expanded in our derivation of the flow equation. However, it quickly becomes apparent that this 
is a systematic error and largely independent of the particular state. We attribute this to the 
scalar part of the higher order 1/j-corrections that were neglected in our derivation. This only 
serves to shift the entire spectrum by some fixed amount relative to the exact results. For a 
more meaningful comparison we eliminate this shift by considering the excited energies relative 
to the ground state, i.e. we consider the excitation energy (or gap) AEn = En — Eq rather than 
En itself. As seen in Figures 14.31 (b) and 14.41 (a), (b) there is indeed a very good correspondence 
between the exact results for AEn and those obtained from -f^e//- From previous studies [HH] it 
is known that the critical point Ac = 0.5 is characterized by the vanishing of the gap between 
the ground state and first excited state. This agrees with the result of the local approximation 
that AEi ~ \/l — 2A. Figure ^31 (c) shows this gap together with those obtained using -f^e//- 

Whereas the local approximation provides reasonable results for the low- lying states, it fails 
to do so for the highly excited states, as is clear from Figure 14.41 (d). For these states the 
assumption that (J^) + j ~ 0{j^) is no longer valid. This means that the local solution within 
a neighbourhood of = — 1 is insufficient since the ffow relevant to these states occur at points 
distant from Xg = —1. 

Figure 14.51 (a)-(d) provides a global view the spectra obtained using the different methods. 
The qualitative properties of the exact result and that of H^jf are clearly very similar. In 
particular we note the level repulsion among the low-lying states, which leads to the vanishing 
of the gap AEi at the critical point. This behaviour is clearly absent in the RWA case, where 
there appears to be no correlation between levels belonging to different sectors. As expected the 
local approximation predicts the correct behaviour for the low-lying states but fails at higher 
energies. In fact, the local approximation predicts a continuous spectrum at A = 0.5. 
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(c) 



(d) 



Figure 4.4: Various excitation energies. In (a), (b) and (c) j equals 100, while in 
(a), (b) and (d) Pmax = 15. The four solid lines in (c) correspond to Pmax = 5, 10, 15, 20 
as they appear from right to left. The bold line corresponds to the exact solution 
for j = oo or, equivalently, the result of the local approximation. In (d) j equals 40. 



The Q-preserving parts of fieff and {Jz)eff may be constructed in the same way as we did for 
the individual submatrices of H^jf. The calculation of the expectation values is now a straight 
forward procedure, the results of which appear in Figure 1131 feVfh). We again find good agree- 
ment with the exact results. To eliminate any constant shift owing to higher order scalar correc- 
tions we consider the expectation values relative to the ground state, i.e. A {Jz)n — {Jz)n ~ {Jz)o 
where (J^)^ = (£'„ | J^l En). 
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(g) (h) 

Figure 4.5: (a)-(d): The energies of the first 55 eigenstates of the Dicke Hamiltonian 
obtained using different methods. For all of these j = 100 and Pmax = 15. (e)-(h): 
Expectation values of relative to the ground state, j = 100 and Pmax = 18 through- 
out. The energies of the first 55 eigenstates of the Dicke Hamiltonian obtained using 
different methods. For all of these j = 100 and Pmax = 

15. 
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Figure 4.6: (a) The function nj, '"^^(x) for q = Q/j = 2 and A = 2.2 in the ^ — > oo 
Hmit. (b) The ground state energy together with the exact result for j = 200. The 
expectation value of (c) and (d) h as functions of the coupling strength. 



In the final part of this section we present the results obtained by diagonalizing the Q-sectors 
through a second application of the flow equation. All the results obtained thus far can be 
reproduced using this method, and we will not restate them here. Instead we focus on the RWA 
case, particularly with respect to the phase structure. The RWA Hamiltonian 



Hrwa = Jz + n + 




exhibits a phase transition at A = 1 .38_ , in contrast to critical value of Ac = 0.5 for the full 
Dicke Hamiltonian. The two phases are distinguished by the order parameter {Eq \ h\ Eq) / j that 
becomes non-zero only in the second phase. This signals a macroscopic occupation of the bosonic 
mode which is responsible for the phenomenon of super-radiance. We proceed by solving the flow 
equation (|4.48|) for a range of A and q values. In the £ — > oo limit the off-diagonal part ni{x,£) 
vanishes and we obtain a set of functions n'^'''\x) = no(x,oo). Figure UHl (a) shows a typical 
example for q = Q/j = 2a.tX = 2.2, together with the exact result for j = 200. Due to the 
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ordering of the eigenvalues at £ = oo we expect each Uq ' (x) to be an increasing function with 
jn'^''^\—l) corresponding to the lowest eigenvalue. The overall ground state energy is given by 

^o=imin{4^''^Vl)}- (4.50) 

The result of this calculation appears in Figure EEl (b). Once it is known to which sector the 
ground state belongs for a certain A we solve the flow equation for the observables h and 
within this sector. Figures 14.61 (c) and (d) show the ground state expectation values calculated 
in this manner. 



Conclusion and outlook 



We hope that our presentation has convinced the reader that flow equations obtained from 
continuous unitary transformations present a versatile and potentially very powerful technique 
for the treatment of interacting quantum systems. Our results for the Lipkin and Dicke models 
have clearly illustrated the myriad of information yielded by a non-perturbative solution of the 
flow equations. This approach is not without its difficulties however. As discussed in Section 
12.21 the construction of smooth initial conditions, although relatively simple for the models 
considered here, are generally non-trivial and in some cases possibly on par in difficulty with 
diagonalizing the Hamiltonian itself. Although the introduction of additional variables may 
aid in this construction, the resulting high dimensional PDE presents its own challenges, as 
was encountered in the Dicke model. Further studies into these issues are needed if the flow 
equations are to become a truly general and robust framework for non-perturbative calculations. 
It also seems inevitable that the treatment of realistic models would require further model-specific 
approximations which were largely absent in our approach. One possibility is the generalization 
of the local approximation scheme of Sections 13.71 and 14.4.21 The flow equation would be solved 
locally around an appropriate point, possibly determined by a variational method, and using a 
specialized generator. This may allow for a treatment of the system by considering fluctuations 
of the degrees of freedom around their classical values. 

The flow equation is one example of a non-linear operator equation of which there are nu- 
merous others appearing in virtually all branches of physics. Although our focus has been on a 
small subclass of these, the solution methods developed here are much more general in nature 
and may in future find application in other fields as well. 
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APPENDIX A 
Representations by irreducible sets 



In Section 12.11 we made use of the fact that any flowing operator may be expressed in terms of 
operators coming from an irreducible set which contains the identity. Now we show why this is 
the case, using a result by von Neumann. First we establish some notation. Let TC he a finite 
dimensional Hilbert space and denote by B{7i) the set of all matrix operators acting on it. The 
commutant Ai' of a set M C B{TC) is defined as the set of operators which commute with all 
the elements of A4. The double commutant of A4 is defined by A4" = {Ai'}'. The following 
theorem provides the key: 

The Double Commutant Theorem: If A is a subalgebra of B{7i) which contains the 
identity and is closed under hermitian conjugation then A = A" . 

Now suppose S C B{7i) is an irreducible set of hermitian operators, one of which is the 
identity, and A the subalgebra of BiTL) spanned by all the products of operators from S. By 
Schur's lemma S' = A! = \ ^ C} and thus A" = B{7i). It follows from the theorem above 
that A = B{7i), i.e. any operator acting on B{7i) may be written in terms of operators coming 
from S. 
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APPENDIX B 

Calculating expectation values with respect to coherent states 



Define 1^) = (1 + 2:2;*) ^\z) where I2) = exp(2j+) \ —j) is tlie unnormalized coherent state. These 
states are known to provide a basis for the Hilbert state through the resolution of unity 

2j + l f dz 

where the integral ranges over the entire complex plane. Furthermore, for any 2 7^ the state \z) 
has non-zero overlap with every state in the Hilbert space. The action of the su(2) generators 
on these states may be expressed in terms of differential operators with respect to 2 and 2*: 

Jz\z) = (-j + 2^^ \z) = X (2, d,) \z) (B.2) 

J+|2) = ^|2) = J+(2,5,)|2) (B.3) 

J-\z) = {^zj -z^^^\z) = L (2,5,) 1 2). (B.4) 

We will always consider 2 and 2* to be distinct variables, and that {z\ is a function of 2* only. 
When calculating expectation values of the double commutators in Section 12.11 one encounters 
expressions of the form 

{z\H^,H^,...Hr.^z) = iAHn.H .Hr^M (3^5) 

\Z\Z) 

where rij = or 1. These can easily be computed by replacing each operator by its differential 
representation to obtain 

{z\Hn^Hn^...HnJz) = -^Hn^ {z,d^) . . . Hn^ {z,dz) {z\z), (B.6) 

\Z\z) 

where the orders of the operators have been reversed. As a result these expectation values can be 
expressed as some rational function of 2 = rexp(i^) and 2* = rexp(— i^). Applying this method 
to (Hq) and (^^i), as defined in Section ITTl we obtain 

,v ./r-2-l\ ,,,,,,, (2i - l)r2 cos(20) ,^ ^, 

(^l^ok) = J [^^) and {z\H,\z) = ^^^^2)2 • (^-7) 
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B. Calculating expectation values with respect to coherent states 
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Using these equations to express r and 6 in terms of {Hi) we may consider expectation values of 
the form of equation (|B.5|) as functions of these averages. For example 

{z\[[AHo, AHi],{AHof]\z) = 8jxy, (B.8) 

where (Hq) = jx and (Hi) = jy. 



APPENDIX C 
Scaling behaviour of fluctuations 



We wish to show that 

(z|Ai?„jAi7„2...Ai7„Jz) ~iL|J for V t G N where G {0, 1} . (C.l) 

From dh'ect calculation this is found to hold for t = 2, and we employ induction to obtain the 
general result. Assuming that this holds for all products of /c < t fluctuations, the induction step 
consists of adding either AHq or AHi to a general product M = AHn^AHn^-'-AHn^ and then 
proving the result for {z\MAHi \ z). 

First we add an extra AHq. Using the results from the previous section we may write 

{z\MAHo\z) = -j{z\M\z) + {z\Mzd,\z) - {Ho){z\M\z) (C.2) 
= -2jzz*{l + zz*)-\z\M\z) - z{z\{d^M)\z) 

+zd,{z\M\z) (C.3) 

where ()B.7|) and {z\Mzdz\z) = z{dz{z\M\z) — {z\{dzM)\z)) were used. It should be remembered 
that M contains a z dependency through the average appearing in each AHi. Dividing by {z\z) 
leads to 

{z\MAHo\z) = -2jzz*{l + zz*)-^{z\M\z) - z{z\(dzM)\z) 



The last term can be rewritten using 



jdz{z\M\z) = dz{z\M\z) + -M— (z|M|z), (C.5) 



(l + zz*)2i ' ' ' ' ' ' {1 + zz* 

which is just the product rule, to obtain 

{z\MAHo\z) = zdz{z\M\z) - z{z\{dzM)\z) . (C.6) 

Note the important cancellation of terms proportional to j{z\M\z). From the induction hypoth- 
esis, and the fact that {z\M\z) is polynomial in j, the first term in ()C.6|) will be of order j'-a-'. 
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The second term becomes 

t 

zY,{z\^Hn,...AHn,^,AHn,+,...AHM{d,{Hn,)). (C.7) 

i=l 

Taking into account the expressions for {Hq) and {Hi), we conclude that each term in equation 
()(y.7|) will have order j^^L— J ^ which reduces to in both the cases where t is odd and even. 

Thus, for any product AI of t fluctuations we arrive at 

{z\MAHo\z) j^'-^K (C.8) 

which concludes the induction step. 

Exactly the same procedure is followed when adding a AHi , although more algebra is required 
as Hi(z,dz) now contains second order derivatives to z. The final result remains unchanged: 

(z|mAi7i|z) (C.9) 

In this case these results make exact the general notion that relative fluctuations scale like 
powers of one over the system size. Finally we mention that when calculating the expectation 
values of a double commutator there is often a cancellation of leading order terms. This can be 
seen in equation (|B.8() . where the sum of terms of order j'^ turns out to be of order j. 



APPENDIX D 
Decomposing operators in the Dicke model 



We define two classes of operators: Cn,m = [J+T i^T + (J-T {b^)"" and r„,„ = (J+)" (^t)"" + 
(J_)" (6)™. Note that Cn.m changes Q hy \n — m\ while Tn^m does so by n + m. We wish to 
find an algorithm for writing a given Cn,m in terms of Tq = J^h + J_6^ and T.^.'s for which 
\n — m\ = n' + m! . Central to this procedure are the relations 

TqCn-i,m.-i - (j^ " Jf) nCn-2,m-2 when m,n>2 

_ I TgTo^rn-i - nTi^rn-2 when m > 2, n = 1 

TqT^-1,0 - {f - Jf) Tn-2,1 when m = 1, n > 2 

Tq when m, n = 1 

which follows directly from the definitions and basic properties of spin and boson operators. 
Given a Cn.m we apply these relation repeatedly, until finally only T.^.'s and Tq appear. Note 
that the former can only appear linearly. Furthermore, any Cn',m' appearing in intermittent 
steps has \n' — m'\ equal to \n — m|, and any Tn'^m' has n' + m! = \n — m\. 
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APPENDIX E 
Dicke model flow coefficients 

The non-zero entries of the A^^^ matrix. For clarity the subscript indices have been raised. 



A 
A 
A 
A 

A' 
A 
A 
A 
A 
A 
A' 
A' 
A 
A 
A 
A 
A' 
A 
A 
A 
A' 
A 
A 



(0) 


(0,Xs) 


{^,Xg) 




-Zlxg (-1 + X5^) 


(0) 


{0,Xb) 






— 16xft 


(0) 


{0,Xq) 


f r\ \ 

[^,Xg) 






(0) 


{0,Xq) 


i f\ \ 

(0,X6) 




-8Xq 


(0) 


{0,Xq) 


(0,Xg) 




4 (-l + 2xf,Xs + x/) 


(0) 


(1) ^JsJ 


{^,Xg) 




-ml' 11 2\ 2 

-16(-l-|-x/) 


(0) 




\^,Xq) 




4 (-1 -l-x/) Xq 


(0) 








-8 (2X6 (-1 +iCs ) -\-Xq ) 


(0) 




(0,Xg) 




-4 Xb Xq 


(0) 




(0,Xs) 




8 (-1 +X/) Xq 


(0) 


i,!) Xq) 


(0, Xq) 




8x6 (-l + a;s ) 


(0) 


(2, Xgj 


(0, Xb) 




-8 (2x6 (-1 + 2;/) + V) 


(0) 


[2, Xg] 


(0, Xg) 




4 (-1 -|-X5^) Xq 


(0) 


(2,2:5) 


(0,Xf,) 




— 16X6^ 


(0) 


{'^,Xb) 


(0, Xg) 




4 X6 Xq 


(0) 


(2, Xq^ 


(0,X6) 




-8X6 Xq 


(0) 


[2, Xq) 


(0, Xq) 




8x6 (-1 + Xs^) 


(0) 


(3, Xg) 


(0,Xs) 




-16Xs (-1 + X^^) Xq 


(0) 




(0,X5) 




8 (-l + X^^) Xq 


(0) 


[3, Xg) 


(0, Xq) 




8 ^2x6Xs (-1 Xs^) - (-1 Xs^)^ XsX 


(0) 


(3, Xg ) 


(l,a^&) 




-4Xq (4X6 (-1 + X^) Xg^) 


(0) 


[3, Xg) 






4 (-1 + x^^) (4x6 (-1 + Xg^) + x^) 


(0) 


(3, Xg 1 


%Xb) 




4xg (4x6 (-1 + X5^) +Xg2) 
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A(o) 


[3, Xs) 




= -4 (-1 + xs^) (4x6 (-1 + Xs'^) + X 


A(o) 


(3, Xb) 


(0, Xs) 


16 Xf) X g Xq 


A(o) 


(3, Xb) 


(0, Xfe) 


= -SXbXq 


A(o) 


(3, Xb) 


(0,Xg) 


= -4 (2xb (-l-2x5X^ + x^2) +Xg2) 


A(o) 


(3, Xb) 


ihXs) 


= 4Xg (4X6 (-1 +X,2) +Xq^) 


A(o) 


(3,Xb) 


(l,a;q) 


= 4X6 (4X6 (-1 + Xs'^) + Xq^) 


A(o) 


(3,2:6) 


(2,x,) 


= -4X5 (4X6 (-1 + Xs'^) + Xq^) 


A(o) 


{3,Xb) 


(2,Xg) 


= -4X6 (4X6 (-1 + x/) + Xq2) 


A(o) 


y3, Xq) 


(0,Xs) 


= -32x6Xs (-1 + Xs^) 


A(o) 


(3, Xq J 


(0,X6) 


= 16x6 (-1 + x/) 


A(0) 


1,3, Xq) 


(0,Xq) 


= 4 (-1 + 2x6Xs + Xs^) Xq 


A(o) 


1,3, Xq) 


(l,a;s) 


= -4 (-1 + X,2) (4.Tfc (-1 + .T,2) + X 


A(o) 


(3, Xq ^ 


(l,a;b) 


= -4X6 (4X6 (-1 + + Xq^) 


A(o) 


1,3, Xq) 


(2,x,) 


= 4 (-1 +Xs2) (4X6 (-1 +X,2) +Xq2 


A(o) 


(3, Xq y 


(2,a;6) 


= 4x6 (4x6 (-1 + a;/) + Xq^) 



The non-zero entries of the A^^^^ matrix. 



A(i) 


(l,0)(0,x,) 


= -4 


A(i) 


(3,0)(0,Xq) 


= -2 


A(i) 


(3,0)(l,x,) 


= 2Xq 


A{i) 


(3,0)(l,Xq) 


= 2X6 


Ad) 


(3,0)(2,x,) 


= -2Xq 


A(i) 


(3,0)(2,Xq) 


= -2x6 


A(i) 


(3,x,)(l,0) 


= 2 Xq 


A{i) 


(3,x,)(2,0) 


= 2Xq 


Ad) 


(3,Xq)(l,0) 


= -2x6 


Ad) 


(3,Xq)(2,0) 


= 2x6 
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The non-zero entries of the A^^^ matrix. 



A(2) 




= -4 




A(2) 


V ' / V ' <ll 


= 4x= 




A(2) 


f3 OVl jjf,) 


= —2 T„ 




A(2) 




= 2 (—1 + X 


^) 


A (2) 




— ^a^g 




A(2) 


(3,0)(2,xg) 


= 2-2xs'^ 




A(2) 


(3,X6)(1,0) 


= 2xg 




A(2) 


(3,Xb)(2,0) 


= -2xg 




A(2) 


(3,x,)(l,0) 


= 2-2x/ 




A(2) 


(3,xg)(2,0) 


= 2 (-1 + X. 





The non-zero entries of the A*^^) matrix. 



A(3) 


(l,0)(0,x,) 


= -2 




A(3) 


(2,0)(0,xg) 


= -2 




A(3) 


(3,0)(0,x,) 


= 8x5 




A(3) 


(3,0)(0,X6) 


= -4 




A(3) 


(3,0)(l,x,) 


= 4(-l + x. 




A(3) 


(3,0)(1,X6) 


= 4xfe 




A(3) 


(3,0)(2,x,) 


= 4-4x,2 




A(3) 


(3,0)(2,X6) 


= -4x6 




A(3) 


(3,X3)(1,0) 


= 4-4x,2 




A(3) 


(3,x,)(2,0) 


= 4(-l + x. 




A(3) 


(3,X6)(1,0) 


= -4x6 




A(3) 


(3,X6)(2,0) 


= 4x6 
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